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(~| Abstract 

-(— > 

^ We consider 2-diniensional Schrodinger operator with the non- 

degenerating magnetic field in the domain with the boundary and 
^NJ under certain non-degeneracy assumptions we derive spectral asymp- 

^ totics with the remainder estimate better than 0{h~^), up to 0{fi^^ h^^) 

—L. and the principal part x /j~ where /? <C 1 is Planck constant and 

Cvl /U ^ 1 is the intensity of the magnetic field; fih < 1. 

^^ We also consider generalized Schrodinger-Pauli operator in the 

^O same framework albeit with fih > 1 and derive spectral asymptotics 

with the remainder estimate up to 0(1) and with the principal part 

X Hh~^, or, under certain special circumstances with the principal 

11 
part X //2/j 2 



cd 0.0 Introduction 



Our goal is to derive spectral asymptotics of 2-dimensional Schrodinger 
operator 

(0.1) A = J2Pjg^''Pk+V, with Pj = hDj - fiVj 

j.k 

near the boundary where 



(0.2) g-''^ = g'^J , Vj and V are real-valued functions, h G (0, 1], /i G [1, oo). 

and in 6(0, 1) C X the following conditions are fulfilled: 

(0.3)i_3 |DV1 < c, \D'^Vj\<c, |D"\/|<c ^a:\a\<K, 

(0.4) eo < ^g^'^VjVk ■ \v\'^ <c Vr/ G M^ \ Vx G 6(0, 1) 

J.k 

So we basically want to generalize results of Chapter 13 [lvr2]^) d = 2 to 
the case of supp-^ n dX ^ 0. We assume that condition (2.60) is fulfilled: 
Fxl. 

However it is not a simple generalization as propagation near the bound- 
ary is completely different from one inside of the domain. While classical 
dynamics inside is a normal speed cyclotron movement combined with a slow 
(with the speed 0(/i~-^)) magnetic drift, it is not the case near boundary: 
when cyclotron hits the boundary it refiects from it and we arrive to a 
normal speed (with the speed 0(1)) hop movement along the boundary. 

The really difficult part is that hop movement is not separated from 
cyclotron plus magnetic drift movement: first, as we move away from 
the boundary the former is replaced by the latter; second, during some 
hop the hop movement can be torn away from the boundary and become 
cyclotron plus magnetic drift movement and v. v.: cyclotron plus magnetic 
drift movement can collide with the boundary and become hop movement. 

The main goal is to investigate the generic case as VVF^-^ 7^ and 
Vex\/F-i = ^ V|x^/^"^ 7^ on dX. 

Plan of the article 

Section 1 is preliminary: first, we consider a classical dynamics, described 
above, in details. Then we consider a model operator in the half-plane and 
derive precise formula for it. 

In section 2 we consider a weak magnetic field case ///? ^ 1 when classical 
dynamics defines everything. Then hop movement breaks periodic cyclotron 



^) This article is a rather small part of the huge project to write a book and is just 
Chapter 15 consisting entirely of newly researched results. Chapter 13 corresponds to 
Chapter 6 of its predecessor V. Ivrii [Ivrl] . External references by default are to [Ivr2] . 



movement which allows us to prove that the contribution of this zone to the 
remainder is 0{fi~^h~^) where factor fi~^ is the width of the boundary zone 
and 1 is the time for which we typically follow classical dynamics. Recall 
that the contribution of inner zone to the remainder (under appropriate 
non-degeneracy conditions) also is 0{fi^^h^^) albeit there factor fi^^ comes 
from time x /i for which we typically follow a classical dynamics. Sure, there 
is a transitional zone between boundary and inner zones but as magnetic 
field is not strong, it is very thin. 

In section 3 we study a strong magnetic field case. In subsections 3.1-3.3 
we establish Tauberian remainder estimates under different non-degeneracy 
assumptions. 

In subsection 3.4 we study propagation of singularities in the transitional 
zone and find that under Dirichlet boundary condition we are able to prove 
better results than under Neumann boundary condition. This difference 
is not technical as it is the first manifestation of the fact that as magnetic 
field grows stronger the classical dynamics loses its value as predictor of 
the propagation and in the case of Neumann boundary condition some 
singularities propagate along the boundary in the direction opposite to hops; 
this is not related to a magnetic drift but rather to a different behavior of 
the eigenvalues of the model operator. 

In subsection 3.5 in the case of the strong magnetic field we pass from 
estimates to calculations and derive our final results. 

In section 4 we consider the cases of super strong magnetic field (/i/j x 1 
and fih ^ 1 respectively; in the latter case we need to take potential 
—lnhF + V rather than V to prevent pushing bottom of the spectrum 
to high. Here 3 > 1 under Dirichlet boundary condition but could be 
smaller under Neumann boundary condition; then spectral asymptotics are 
concentrated near boundary). As we know from Chapter 13 [Ivr2] cyclotrons 
are no more observable due to uncertainty principle but magnetic drift 
preserves its sense. The same remains true near boundary: we do not see 
hops but we observe propagation along the boundary and it can be torn 
away from the boundary and became magnetic drift and v. v. 

Finally, in section 5 we consider generalizations: first, we get rid off 
condition V x 1 and then we get rid off condition F x 1. In the latter case 
as we know from Chapter 14 [Ivr2] our remainder estimate must deteriorate 
to 0(/i-i/,-i). 

Appendix 6 is devoted to auxiliary 1-dimensional harmonic oscillator 
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—d'^ + x^ considered on (— oo.r^] with Dirichlet or Neumann boundary 
condition as x = rj. We study its eigenvalues Xo.niv) and X^.niv) reproducing 
results of B. Helffer, C Bolley and M. Dauge. These properties are crucial 
for our formulae and for analysis of section 4. 

1 Preliminary discussion 

1.1 Inner and boundary zones 

Let us consider 2-dimensional Magnetic Schrodinger operator near the 
boundary. Let as usual 7(x) = |dist(x, 9X) and consider disk (i.e. 2- 
dimensional ball) B[y, 7(y)). Scaling it x i— > (x — y)'y~^ to the unit disk we 
get 

(1.1) /? h^ /?new = /'7~^ /i tH- /Unew = /i7. '^ ^ ^^new = «^7. 

i^ ^ i^new = 1>1 

where recall v and D are introduced in (13.3.11) and 13.2.77 [Ivr2]. We 
assume that originally v = D = 1. Then according to Chapter 13 [Ivr2] 
under reasonable conditions contribution of B[y, 7(y)) to the remainder 
R^^ does not exceed 

as 'J > fih and dividing by 7^ and integrating we conclude that the contri- 
bution of the inner zone 

(1.3) Xinn:={x,7(x)>7:=Co/i-^} 

to R^^ does not exceed 



;i.4) Cfi-^h-^ f -f ^ dx X Cfi-^h-^ log/i 



as /i < /? 2 . 

Remark 1.1. We must take Cq large enough in (1.2) to ensure that /Xnew > A^o 
where /iq is large enough constant. 
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One can get rid off logarithmic factor in (1.4) using Seeley' approadi 
(R. Seeley, [SI, S2]; see also section 7.4) [Ivr2]. One also should consider 
case fi > h~2 but then zone {fih > 7 > Cofi~^} appears and its contribution 
to remainder is estimated differently, as 

(1.5) C f -f{x)-^T{x)-^dx. 

We will postpone this until real proofs. 
In the boundary zone 

(1.6) Abound :={x,7(x)< 27} 

one can apply standard Weyl remainder estimate after rescaling. Then 
contribution of B{y, 7) to R^ does not exceed Ch~^^ = Ch~^^ and the total 
contribution of Xbound to the remainder does not exceed Ch~^^ x ^~'^ = Ch^^. 
This is much worse than the contribution of the inner zone we derived 
and needs to be fixed. To do it we need to extend dynamics from time 
T X 7 X fi~^ ( T X 1 after rescaling) to 7" x 1 [T ^ ^ after rescaling) or 
even more in the strong, very strong and superstrong magnetic field cases. 

1.2 Classical dynamics near the boundary 

To understand the role of the boundary, consider classical dynamics. Let 
us consider first a half-plane X = M^ = {->^, -'^i > 0}, g'-''' = Sj^, F = 1 and 
V = 0. Let us write the model operator in the form 

(1.7) A:= h^Df + {hD2~fixif, 

as F12 := ^2^/1 — 9i \/2 = —1 according to (13.1.7) [lvr2]. 
Then we have a Hamiltonian circular trajectory 

(1.8) xi = fi^^^2 + 3n^^ cos2fxt, -x'2 = ^2 — 3fi^^ sin 2fit, 

^1 = -asin2/^t, 6=6 

where a = t^ and r is an energy level. So we got circular counter-clock- wise 
trajectories of the radius rs centered at x with xi = /i~"^6 and depending 
on 6 these trajectories behave differently: 
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Figure 1: Different classical trajectories in half-plane for model operator 



(a) As ^2 > T2 trajectory does not intersect dX or just touches it and 
remains circular. 

(b) As .^2 < ^2 trajectory reflects from dX and we get a "hop" -movement: 

(c) As .^2 \ — T2 trajectory stays closer and closer to dX and becomes a kind 
of gliding ray in the limit. And we are interested only at zone {^2 > — ra}. 

So, trajectories described in (b) are not periodic even for the model 
operator. 





(a) 6 > (b) 6 < 

Figure 2: Calculating the length of the hop; p = l — r], p' = l + ri 



One can calculate easily that 

(1.9) The length of the hop (i.e. the distance between hop's start and 
end points) is 2afi~^{l — rf^i with 77 = ^2/^ while the length of the arc is 
2afi~^{TT — arccosT]) and thus the time of the hop is /^"■'"(Tr — arccosrj). 
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Therefore 

1.10) As \r]\ < 1 the average hop-speed along X2 is ^ 

(1-772)5 



-23v[ri) with 



(1.11 



v{r]) 



71 — arccos[ri) 



This is a plot of v[ri). One can see easily that v[ri) is 
defined on ( — 1,1) where it decays from to 1; 77 = 1 is a 
threshold between circular and hop- movement and rj = —1 
corresponds to gliding rays. 

As we cannot get remainder estimate better than 0{^h^^) for model 
operator we need to consider a perturbation by a potential: 



:i.i2) 



A := ItDl + (/7D2 - ^ix^Y + V{x); 



a classical dynamics for a general operator (0.1) in dimension d = 2 will be 
not different in our assumptions. 

Then one should use the same classification as before with 



(1.13) 



V 



6 



Wo{xy 



1/1/ = (r - V{x))F{x)-\ 1/1/0=1/1/ 



dX 



where so far F[x) = 1. 

Actually the last statement is not completely true in the transitional 
zone (where now we consider zones in (x, ^)-space) 

(1.14) A;,3ns = {(x,0 : 16 - H/o(x)5| < 2p} 
or equivalently {jr] — 1| < Qp} where so far 

(1.15) P=Cofi-^ 

but later it may be increased due to uncertainty principle. 
So, let us introduce an inner zone 



(1.16) 



'^inn = {(x,0: 6-H/o(x)5>p} 
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and a boundary zone 

(1.17) A'bound = {(x,0 : 6 - l/l/oW^ < -p}. 
There is no need to consider gliding zone 

(1.18) A'giid = {(x,0: 6+W/o(x)^<p} 

separately from A'bound. 

Recall that inside of the domain potential causes magnetic drift 

(l-19)i,2 j^xi = fi-'d,, W, -X2 = -fi-'d,, W. 

Let us analyze what happens near the boundary. Note first that 

(1.20) Billiards do not branch as /i > /ig where po is large enough. 
Really, one can prove easily that 

(1.21) With respect to Hamiltonian trajectories dX is strongly concave in 
the gliding zone and strongly convex in the transitional zone (and domain 
X has opposite property) as /i > /iq. 

Recall that (according to Figure 2) p = 1 — r], p' = 1 + r], r] = .C2l^o~^. Then 
along Hamiltonian trajectories 

dp' d , -ix -i 1 -3 dxj 

l%.,l/l/o"^(l - ^6W/o"'^) mod 0{p-^fl) 

as xi = 0{p^^p') and therefore for one hop 

S-^"'=^oHl-^6^o-'t) 

with Wq, 1/1/0x2. 'C2 calculated in the middle of it 

= Wo.M'{'i- + Vv{v))v{vr' = ^pWo^M' mod 0{p'{p-' + p')) 

as rj = p' — 1 and 1/(77) = 1 ~ |p' mod p'^. 
Therefore 
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(1.22) Along trajectories of the length < 1 p'exp( — ^ ^Vq^) remains constant 
modulo 0{p'{fi''^ + P'))- 

Remark 1.2. The similar statement would be completely wrong for p because 

as p ~ ^ ~ — l%x2l^o ^ ^"^^ therefore as l/l/ox2 < hop-trajectories^-* will 
be torn out of the boundary and begin magnetic drift movement. Meanwhile 
as l/l/ox2 > trajectories drifting in the inner zone may collide with the 
boundary and begin hop- movement ^^ 

In other words hops move away from the boundary (to the boundary) in 
the direction along the boundary, in which Wq decreases (increases). 

Example 1.3. Meanwhile W^^ has more subtle effect. As hop-speed is larger 
than CqP^^ (i.e. in A'bound) magnetic drift with respect to X2 has no qualitative 
effect. However there are no hops in X\„„. Therefore as l/l/x^ x 1 we have 
two rather different cases: 

(i) 1/1/xi > 0. Then according to (1.19)2 magnetic drift is to the left, in the 
same direction as hops. Then all dynamics is to the left^-*. In particular 
as l/l/x2 > hop-trajectories are torn from the boundary and begin drift 
movement (see figure 3(a)) while as W^^ < drift-trajectories collide with 
the boundary and begin hop- movement (see figure 3(b)). 

(ii) H/xi < 0. Then according to (1.19)2 magnetic drift is to the right, in the 
opposite direction to the hops. So direction of dynamics (with respect to 
X2) in A'inn is opposite to the hop-movement. In particular as l/l/x2 < — eo 
hop-trajectories are torn from the boundary and begin drift movement (see 
figure 3(c)) while as W^^ > cq drift-trajectories collide with the boundary 
and begin hop- movement (see figure 3(d)). 

This is consistent with the fact that drift trajectories are level curves of 
1/1/. 

Example 1.4- Assume now that l/l/xj vanishes at some point but 1A^X2X2 7^ 0- 
Then repeating analysis of the previous example we arrive to the following 
four pictures: 

Again this is consistent with the fact that drift trajectories are level 
curves of 1/1/: in cases (a), (d), (b)-(c) point x is a local minimum, maximum 
and minimax respectively. 



^■^ In the positive time direction (and negative X2-dircction) . 
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4- 




(a) l/14i > 0, 144, > (b) 1/14, > 0, 144, < (c) W^, < 0, 144, > (d) I44i < 0, 144, < 

Figure 3: To example 1.3. Bold lines show hop movement and thin lines 
show drift movement which is along level lines of 1/1/. 




(a) W,^ > 0, I44,x. > (b) 144, > 0, 144,x. < 




^^^-^ 



(c) 144, < 0, 144,., > (d) 144, < 0, 144,x. < 

Figure 4: To example 1.4. Bold lines show hop movement and thin lines 
show drift movement which is along level lines of 1/1/. 



However the following observation basically remains true: 

(1.23) The speed of magnetic drift is 0{fi~^) (the typical speed is x fi^^] 
while the speed of hop-movement is 0(1) (and the typical speed is x 1). 



1.3 Spectrum of the model operator 

Consider model operator (1.7) as xi > with the Dirichlet or Neumann 
boundary condition tv|xi=o or C>itv|xi=o- Making /)-Fourier transform we arrive 
to 1-dimensional operator 



^1.24) 



h^Di + {^2-i^x,y 



Chapter 1. Preliminary discussion 11 

which after transformations Xi i— t- fxxi or Xi i— ;■ hr^xi with h := yU^-*-/) becomes 

(1.25) h^Dl + {i2-xif with n = /x/? 
or 

(1.26) L{ri) := h{Dl + (r/ - xi)^) with 7] = h~'^2 

at ]R+ 9 Xi respectively, again with the Dirichlet or Neumann boundary 
condition at xi = 0. As in the previous chapter one needs to distinguish 
between h = fih and fi = ii~^h. 
Obviously 

(1.27) Each of these operators has a simple discrete spectrum, let < 
A*,o(^) < A*, 1(77) < A*2(^) < ■■■ be eigenvalues of operator L^^q) defined by 
(1.26) a,s h= 1 where * means either D or N and as usual D and N denote 
Dirichlet and Neumann respectively. 

Further 

(1.28) Let t'*j(xi, 7]) be real- valued orthonormal eigenfunctions of operator 
/.*(?7) (as h = 1) corresponding to eigenvalues \*j{ri). 

We will analyze them in details later (see Appendix 6), so far let us 
notice only that X*j{ri) — )■ +00 as 77 — )■ —00 and as A*j(?7) are analytic we 
conclude that 

Proposition 1.5. (i) Spectrum of model operator (1.7) in {xi > 0} with 
Dirichlet or Neumann boundary conditions is absolutely continuous and 
occupies [h ■ inf^ K,o{v)' +00); 

(a) Schwartz kernel of its spectral projector is 

(1.29) e{x,y,T,i2,h) = 

j>Q -^ 

Setting X = y, subtracting /)~W^^(r, h) with A/'^^(r, h) defined by 

(13.1.9)2 [Ivr2]: 

(1.30) Ar^^(r, /i/j) = 

#{j G Z+ : (2j + l)^^hF +V<t}- {2n)-'^{fih)F. 
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and integrating with respect to xi we arrive to the following formal (at least 
for now) equality 

(1.31) f (e{x,x,T,i2,h)-h-^Af^^{T,h))dxi = 

(27r)- V/)-' rY^ih^^- hKj{v)Kj{n-'^^i. v) dv- 

e(r-(2j + i)n)) dx, = h-'^f:'Z.,{r,^^h) 

with 

(1-32) ACZnd(^,^): = 

(27r)"^ /" 5^( / e(r - hX,j{v))vlj{x,. r/) dv - B{t - (2j + l)h))h^ dx, 
Jo j^o -^ 

where in the last transition we rescaled xi. 

Recall that Xj = (2j + 1) are eigenvalues of operator (1.26) on the whole 

line M 3 Xi as /i = 1. 

The following properties of X*j{ti) are useful to know 

In virtue of proposition 6.1 r > (2_/ + l)h <^==^ r > \dj{v) as < r^ is 

large enough and thus 

(1.33)d e(r - XdAv)^) = e(r - (2j + 1)^) ^= < 77 > 1 

However r > {2j + l)h <^=^ r > \hj{j]) as < ?] is large enough and in 
this case 

(1.33)^ e(r - XnAv)^) = Q{r - (2j + 1)^+) ^= < 77 > I 

Recall that Q(t) is semi-continuous from the left and therefore 9(0) = 0. 
Meanwhile Q{t+) is semi-continuous from the right and therefore 9(0) = 1. 
Thus for Dirichlet boundary problem formulae (1.31)-(1.32) should be 
retained 

(1.31)o f{eo{x.x.r.h) - /)-W^^(r, n)) dx, = h-'AT^ZUr. ^^h) 
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with 

Jo j^Q V / 

but for Neumann boundary problem they need to be adjusted to 

(1.31), P {e^{x,x,r,h) - h-^M'''^{r+,h)) dx, = h-'^fff^^^„,{r, fih) 

with 

(1-32), <lnd(r,^):= 
(27r)-^y""^(y"e(r-fiANj(ry)X.(xi,r/)c/r7-e(r-(2i+l)n+))n^c/xi. 

Proposition 1.6. (i) Equalities (1.31)p-(1.32)[5 and (1.31)|^-(1.32)|^ hold 
and all integrals converge; 

(a) Alternatively 

(1-34)d KZUr.h) = 

{2n)-' 5^ /(e(r - hXojiv)) - e(r - {2J + l)h)d{v))h^ d^ 

and 

(27r)-^ Y.I[Q{r- hXnAv)) -Q{r- (2j + l)/i+) e(r^)) ^^ c/r^ 
i>o 

Proof. One can prove easily that 

(1.35) |xi - r^l > e(xi + Ir^l) ^ |t;,,y(xi, r^)] < C,-(l + Xi + |^|)-^ 

By definition 



(1.36) 



/■oo 

/ vlj{xi,ri)dxi = l 
Jo 
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and based on these two facts and (1.33) ^ |^ one can prove (i) easily. 

Inserting 7] in the second term in the big parenthesis in (1.32)p ^ through 
factor 1 = / 5(xi — r/) drj, changing order of integration and calculating 
integral with respect to xi we prove (ii). D 

Remark 1.7. From forthcoming weak magnetic field arguments it follows 
that in the weak sense 



(1-37) KToUr.h)^}_^K^,{r)!f'^ 

n>0 

1 

with K*o = =F(4vr)~-'^r_^ as * = D, N respectively. 

2 Weak magnetic field 
2.1 Precanonical form 

We will consider a general magnetic Schrodinger operator (0.1) in X C M^ 

(2.1) ^= Yl PjS^^Pk + ^' with Pj = hDj - fiVj 

satisfying in 6(0, 1) C X assumptions (0.2)-(0.4) but in contrast to (13.1.5) 
[Ivr2] we assume that 

(2.2) X n 6(0, 1) = {xi >0}n 6(0,1). 

Then in contrast to Chapter 13 we cannot unleash a full power of Fourier 
integral operators as we must preserve boundary, but we can assume that xi 
is small, in fact as small as fi~^Pj are: not exceeding C max(/i~-^, /i"?/;?""^). 
Without any loss of the generality one can assume that 

To achieve subsequently (2.3)j^_3 we just reintroduce xi := d\stg/F{x,dX) 
in the given metrics (g-'^^F^^), then reintroduce X2 := a{x) + /3(x2) with an 
appropriate function a(x) and arbitrary function /3(x2) and finally chose 

/3(X2). 
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Remark 2.1. In this construction F is any arbitrary positive function but 
we select it to be a scalar intensity of the magnetic field. We cannot however 
choose c^xig^^lxi^o as this is defined by the curvature of the metrics {g^^F~^). 

Then 

(2.4) Fi2 1x1=0 = -1 

(actually it may be 1 but then we change X2 i— > — X2) and without any loss 
of the generality one can assume that Vi = and V2 = xi + O(x^); we can 
always reach it by a gauge transformation. 
So, changing V by 0{h^) 

(2.5) /l = (/)Di)F(x)(/)Di) + 

(/7D2 - /ixi - iixlb[x))g''\x){hD2 - /ixi - iixlb[x)) + V{x) 

with V(x) satisfying (2.3)3. ^^^ can easily generalize to such operator 
results of the previous section. 

2.2 Propagation of singularities 

2.2.1 No-critical point case 

In the case of a weak magnetic field 

(2.6) /io < /i < h^'^ 

we know that inside of the domain the shift during the first winding is 
microlocally observable under condition (13.3.54) [Ivr2] i.e. as 

(2.7) \VVF-^\ > to. 

We are going to prove that near boundary the same is true as 

(2.8) \VdxVF-^\>to 

where V^x is a derivative along boundary (i.e. 8x2)- 

To prove this assertion one needs just to prove the following statement: 
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Proposition 2.2. (i) The propagation speed with respect to {x,C,2) does not 
exceed C ; 

(a) Under condition (2.8) for operator in the form (2.5) propagation speed 
with respect to ^2 is of magnitude 1 as xi < Cfi~^. Namely, for time t shift 
with respect to ^2 is of magnitude \t\ and has the same sign as tW^^ with 
W = {t- V)F-\ 

Proof. Proof follows the very standard way of the proof propagation of 
singularities as in theorems 2.1.2 and 3.1.2 [Ivr2]. Here we are using an 
auxiliary function (p — ip[x, ^2, t) and not invoking actually reflections from 
the boundary. 

Then proof of statement (i) is then completely trivial. 

To prove (ii) note that modulo 0(xi) 

{a,6} = {F.^2}{3-V) + {\/,6} -2(6 -/xxi -/xx2 6)F/xx2{6,6} 
and the last term is 0{^^^) as \a\ = 0(1); so 

(2.9) {a, 6} = {F, 6}(a - r) + F'{{V - r)F-\ 6} + 0{^^-') 

and as [a — t) is small everything is defined by the second term in the right 
hand expression. 

Further details of this very simple proof are left to the reader. D 

Then we immediately arrive to 

Corollary 2.3. Let U = U{x,y,t) be a Schwartz kernel of e~'^ ^^. Let 
ip{x) = tp' [x2)ip" {fJ^Xi) with fixed functions ip', ip" G ^q°° . 

Finally, let condition (2.8) be fulfilled on supp ■?/''. Then under condition 
(2.6) both 

(2.10) F,^,-i,XT{t)r{U^) 
and 

(2.11) F,^,-^,{xT{t)-XT'{t))^{U^|J) 

are negligible as eifi^^ < T' < T < t2 with arbitrarily small constant e-i and 
small enough constant t2 ^-^ . 



3) Recall that as usual x & '4°°([-l. -5] n [i, 1]), x & '^'^([-l. 1]) and equals 1 at 
'- -] 

2' 2'- 
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Meanwhile rescaling x h-> x^, h ^-^ fih we immediately arrive to 

Proposition 2.4. Let condition (2.6) be fulfilled and let condition (13.2.45) 
[Ivr2] i.e. 

(2.12) V < -eo 

he fulfilled on supp ■?/''. Let ip be specified in corollary 2.3. 
Then as \t\ < e and T = e/i""*" 

(2.13) F,^,-vxr(t)r(ty^) = 

n>0,m>0 n>0,m>0 

where Hnm{T) O'^d n'nmi'^) ^'^^ smooth coefficients. 

Combining proposition 2.4 and corollary 2.3 we immediately conclude 
that 

Corollary 2.5. Under conditions (2.6), (2.8) and (2.12) decomposition 

(2.13) holds with T = €2. 

Therefore applying Tauberian theorem we arrive to estimate (2.14) below 
as ip specified in corollary 2.3. 

Theorem 2.6. Let conditions (0.1)-(0.4), (2.2), (2.6), (2.8) and (2.12) be 
fulfilled on suppV^ where ^p{x) G "^^(R^). The 



len 



(2.14) R^:=| /"(e(x,x,0)- Y. '^nm{0)fi'"h'"^'-')^{x) dx- 

•^^ n>0,m>0 

"^^^ n>0,m>0 

where dSg is a measure on the boundary corresponding to metrics g. 

Proof. So, a,s ip = ilj'{x2)ip"{fiXi) estimate (2.14) is proven; thus contribution 
of zone {x, xi < 2Cofi^^} to the remainder is 0{fi^^ h^^) . This is rather a 
broad zone as Q is arbitrarily large and in zone {x, xi > Cq/x} analysis is 
almost as if there was no boundary. We need to use the time direction in 
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which ^2 increases and therefore trajectories drift inside of X. Still it is not 
quite as without boundary as we are forced to use cut-off functions which 
are scaled with respect to Xi. 

However, this is technicality and we overcome it in the following way: 
first note that such scaling functions are admissible as long as microlocal 
uncertainty principle fi~^ x /i~^ > fi~^h^^^ holds as we reduced our operator 
to microlocal canonical form inside of domain. Here the first factor /i~-^ is 
due to the scale. So, as /x < h^~2 we can appeal to theory of Chapter 13 
[Ivr2]. 

However we can change variables xinew = ^i + bxf so that 

^1 = (l + 0(xi))^inew, 6 - fJ'Xl - fJ'bXi = 6new " /"Xinew + 0{x^)^i 

and then bad term does not appear at all; in old coordinates it would amount 
to replacing ^2 by ^2 + O(x^^i). 

This proves that in the inner zone we can take T x /i in the correct 
time direction and most importantly, we do not scale with respect to X2 
so microlocal uncertainty principle would be 1 x /i~^ > fi~^h^~^ i.e. in our 
frames /i < h^~^. Formula (2.14) is proven. D 

Finally we need to pass to magnetic Weyl formula: 
Theorem 2.7. In frames of theorem 2.6 

(2.15) R^^:=| /"(e(x,x,0)-/)^W^^(x,0,/i/7))7/>(x)c/x- 

JdX 



-u-l 



where A/'J^b^und is introduced by (1.32)^ or (1.32),^ for Dirichlet or Neumann 
boundary condition respectively with h = fihF[x) and r replaced by —V[x). 

Proof. To pass from (2.14) to (2.15) let us recall that coefficients are obtained 
by successive approximation method and notice that as ip is fixed (non- 
scaled) function running successive approximation method to derive (2.14) 
leads to an extra factor h rather than fih if we differentiate ip or S'''^ o^^ ^ o^^ 
differentiate twice Vj; however exactly one such differentiation leads to in 
the final calculations, so we conclude that at least 2 "losing" differentiations 
must be there, extra factor is 0{h^) and the term is 0(1). 
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Thus to derive R^ modulo 0(1) we should not differentiate g-''^, V.ip at 
all and to differentiate Vj only once. However it means exactly considering 
model operator A = Ay ai each point y. 

The microhyperbolicity with respect to X2 implies that with integration 
(with the scale 1) over X2 we can replace TV''^^ by its decomposition with 
respect to powers of. h = fih which justifies the transition from (2.14) to 
(2.15). D 

2.2.2 Analysis in the boundary zone 

In what follows "formula" and "remainder" mean Tauherian formula with 
T = efi~-^ and the corresponding Tauherian remainder R""" until we will pass 
to Weyl and magnetic Weyl formula and remainder. 

To cover the case when condition (2.8) is violated (with an extreme case 
1/1/ := VF^^\qx = const) let us consider first boundary zone A'bound defined 
by (1.17) where we chose a small parameter p later. To do so let us consider 
a stripe 

(2.16) A-bound.p = A'bound H {p < -^ + {VP-^)'^ < 2p}. 

Then the length of the hop is x fi~^p^ with the possible perturbation 0{fi~^) 
due to the magnetic drift (so p ^ p^^ would suffice, but we request p > p~^ 
anyway) . 

Now, weak magnetic field approach would mean that uncertainty principle 
is fulfilled after the first hop: p x p~^p^ > h^^^ or equivalently p > p^h^~^. 
Combining with restriction p > /i"-^ we arrive to 

(2.17) p> p := max(Cop"\p5/?5-^) 

and our goal is to prove that under this the contribution of A'bound to the 
remainder does not exceed Cp~^h~^. 
To achieve this goal let note first that 



(2.18) / dxid^ : da< Cp 

-^A'bound,pn{x2=COnSt} 

which enables us to estimate 

(2.19) \F,^,-i,XT{t)r{u'Qyn < Cp-'h-^p^ 
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as \t\ < e, Q = Q{x2, /JD2) £ S/,,p,^(M) and with the symbol supported in 
'^bound.p n {|x2 - X2I < 7}, p > C max(7, jj,'^), ip = ■?/'"(/iXi) and T = T^ = 

Remark 2.8. Note that factor ii~^ appears here and in ah similar estimates 
because width x p with respect to ^2 matches to the width x pT^p with 
respect to xi on energy levels below c. 

What we need is to investigate propagation until time T = T*{p) and 
prove that 

(2.20) \Ft^hMxM^)-XTM)nuVyn < C'h' 
with 

(2.21) j T*-\p)dp<C 

then as the main part is given by Tauberian formula with T = T^, remainder 
does not exceed 

Cp-^h~^ f T*~\p)dpdx2 < Cp-^h'\ 

In what follows T*[p) = p^^^ with arbitrarily small exponent 5' > 0. 
To prove (2.20) with T*{p) = p^'^' note first that 

(2.22) If |(9x2 VF"-*^! < cp^ then propagation in any time direction remains 
in the zone 

(2.23) {^-p < -6 + {VF-^)\ < p'] n {|X2 - X2I < 27} 

for time T = T*{p) 
and 

(2.24) If l^xjVF^-^l > ^cp^ then propagation in an appropriate^"* time 
direction remains in (2.23) for time T = T*[p); 



^* In which {VF ■'■)(0,X2) decays; as rotation is countcr-clock-wise we take ±t > as 
±(\/F-i)(0,x2) >0. 
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in both cases 7 = Cfi^^ T*{p') < ep for sure where in the former case p' = 3p 

and in the latter case p' = p^^^ . 

After assertions (2.22) and (2.24) are proven, we need to prove that 

singularities really propagate. Figure 2(a) shows that as VF~^ = 1 and the 

movement is strictly circular the length of the hop is exactly 2/i^"'"(l — (1 — 

1 J 

p)^) ^ while in general case it will be of the same magnitude p~^p2 as p <^ 1. 

But obviously in the most critical zone movement in X2 is not monotone. 

However for sure X2 increases with each hop as we bound ourselves with 

xi < ^'f^~^P with small enough constant e'. 

The main problem are trajectories which are almost tangent to the 

boundary. There are two kind of them: with p' ^ 1 and with p ^ 1. 

Trajectories of of the first kind are not actually difficult (see proof of 

proposition 2.11). So we start from trajectories of the second kind. 

Proposition 2.9. Consider zone Z := {,^2 > ~^}- Consider point (3c, ^2) G 
Z and its (7, a)-vicinity Q with 7 = (71, 72) and^^ 



(2.25) 7i = a = m-\ as xi > m-\ 

q 1 c 1 

(2.26) 7i = (T = (/i7"2^2 ^^ ^g xi X 7 < /i2 

(2.27) 72 = h^~^ 



and operator Q = Q{x, hD2) G S/,,^-^,^2,o- with symbol supported in Q. 

Further, consider point {y,fi2) G Z and its {-y, a) -vicinity Q.' with 7 
reintroduced for this point according to (2.25) -(2.27) and also operator 
Q' = Q'{x, hD2) G Sh,'yi,'y2,<T with symbol supported in Q! . 

Let T X p^-*- and 

(2.28) |xi - 7117,-^ + |X2 - y2|72~' + 16 - ^2!^"' > e 

Vx G M/t,^(Q), y eO! yt e suppxr 

where ^t.r is a Hamiltonian flow with reflections and ^1 is defined so that 
a(x, ^) = T. Then 

(2.29) F,^,-^,XT{t)Q'U'Q' 

is negligible. 



^' After rescaling x i-^ jix. 
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Proof. As we assume ^2 > —^0 we do not need to consider "short and low" 
hops as on figure 2(b) and the time of the hop is x fi^^. 

Obviously as /x < h~^ one can apply results of sections 2.4 and 3.5 [Ivr2] 
and justify our final conclusion. 

Let us consider larger /i. To do so we need to understand how small 
vicinity of Hamiltonian billiard flow with reflections we must take to contain 
propagation. 

To do so consider propagation in the different zones, step by step, and 
we do consider (7,71,72 not necessarily defined by (2.25)-(2.27). 

(a) First of all, as xi > 7 > h2~^ one can take any 7^ > 7 and pj > 7. So, as 
xi > h2~^ we can take 71 = h2~^ and other scales different but larger. Thus, 

(2.30) Let xi > 7 > h2~^ and let ■yj > 7, a > 7. Let t be such that 
M/t'(supp (?) does not intersect 7- vicinity of dX as < ±t' < ibt (if it > 0). 

Then (2.29) is negligible, provided v|/t,^(supp (?) does not intersect (7, a)- 
vicinity of supp (?'. 

We refer to this as inner propagation (see figure 5(a)). 

(b) Consider now zone Cp < xi x 7 < h2~^. Let us scale xi 1— )■ xi7~"'^, 

1 3 

X2 ^-)■ X27~2 and H \-^ h' := hj~2. As after original rescaling x 1— )■ /ix we 
had that derivatives of all coefficients were less than Cfi~^ we conclude that 
after this rescaling they are less then Cfi~^'j2 and after division by 7 they 
are still bounded. 

More precisely: we recalled that "up to perturbation" operator was 

(2.31) A = {HDif - {hD2 - xif - 1 

with 1^2 — 1| < p; so in the zone in question rescaling as described is justified. 
And therefore respectively we need to scale ^1 i-> ^i7~2 and (,2 — I ^-^ ^2l~^- 
Now we have rather regular situation and can take h'2~^ vicinity. In 
other words, we can take 

(2.32) 7i,a > h''2~^-f = (^7^5)^"%, 

(2.33) 72 > ^'5-^7^ = {h^-ly'\-2 



and we can replace these right hand expressions by h2~^'^^ and h2^^j^^ 
respectively and the latter does not exceed its value as 7 = p and it is h^^^. 
This is squeezed inner propagation (see figure 5(b)). 
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(a) Inner propagation 



(b) Squeezed inner propagation 




(c) Squeezed propagation with reflection 

Figure 5: Inner, squeezed inner (and similar mirror symmetric), squeezed 
with reflection propagations 



Furthermore, we can apply long-range propagation of section 2.4 [Ivr2] 
and replace condition xi x 7 by a weaker condition 7^+*^ < xi < 7^"^ with 
6' > small enough. But then we need no more than C{S') jumps to reach 
from 7 X ^12^'' to 7 X h^^~^ or inversely. 

(c) Finally as xi < Cp we can apply the same scaling with 7 = p and note 
that after rescaling trajectories meet the boundary under angle disjoint from 
0. So we have a standard reflection situation. Scaling back we arrive to the 
same conclusion as in (b). This is squeezed propagation with reflection (see 
figure 5(c)). D 

Repeating N :< fiT times we arrive to 

Corollary 2.10. Consider fi^^ < T < ep. Then conclusion of proposi- 
tion 2.9 remains true if we redefine a := ComN, ■Jj := 7;,oidA/ with N ^ pT 
number of rotations. 

Now we can prove our main claim 
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Proposition 2.11. Estimate (2.20) holds with T^ = tjj, ^ and T*{p) = 

Proof. We consider a pseudo-differential partition of unity. Let Q' = 
Q'{x2,hD2) G S/,,p,^ witfi symbol supported in e(7, p)-vicinity of (x2,^2) 

and 7 = p'^p^ satisfy (2.17). 

(i) Consider first ^2 < — eo- 
Note first tliat 

(2.34) As ,^2 < — e botfi operator and boundary value problem are ^2- 
microhyperbolic as Xi > 0. 

Really, ih^^[A, X2] = hD2 — pxi < — e. 

Then due to results of chapter 3 [Ivr2] U ^Qy is negligible as ±t > and 
±(x2 — X2 + 67) < ±eof) which implies (2.20). 

(ii) Consider now ^2 > — eo- Note that M/t after exactly one turn moves X2 
to the left by at least epz; we assume that p > Cop~^ to counter possible 
perturbations. Further, as p > {ph)3~^ it exceeds 72 expansion due to 
uncertainty principle. This justifies conclusion of the proposition as T* = eop 
as propagation speed with respect to ^2, ^2 does not exceed Cq. 

(in) To increase 7"*(p) notice that we can take T*{p) = p^~^ unless i > p^ 
where 

(2.35) i:=eo\V9xVF-'\ 

as the speed of propagation with respect to p does not exceed i and dynamics 
remains in the same {i, p)-element with respect to (x2, ^2)- 

Meanwhile as i > p^ we can take time direction in which p increases 
and then take 7"*(p) = p^~^ anyway. 

Uncertainty principle pi > h^~^ is obviously satisfied in both cases. D 

Remark 2.12. Surely one can improve arguments of (iii) and we will do it 
studying strong magnetic field. However at this moment this leads to no 
improvements in the remainder estimate. 

Then immediately we arrive to 
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Corollary 2.13. Contribution of zone A'bound defined by (1.17), (2.17) to 
the remainder does not exceed Cfi^^h^^. 



Proof. Contribution of A'bound.p to the remainder does not exceed 
(2.36) Cfi'^ph'^ X T*--^ X Cfi~H'^p^' 



and summation over p > p (i.e. integration with respect to p ^dp) results 
in 0{p~^h-^). U 

2.2.3 Analysis in the transitional zone 

Consider now transitional zone defined by (1.14). Obviously we get a rather 
rough estimate 

Proposition 2.14. Contribution of the transitional zone (1.14) to the re- 
mainder does not exceed 

(2.37) Cph~^ X Cp-^-^ + C{ph)h-^-^. 

In particular, as p < /?~5+<^ it does not exceed Cp^^h"^. 

Now we want to improve this estimate under some non- degeneracy 
condition invoking VF^^\gx- 

Let us introduce ^-admissible partition with i defined by (2.35). 

Proposition 2.15. On (.-element with 

(2.38) £>£:=max(C/i~\/i5/)5-'5) 

expressions (2.10) and (2.11) are negligible as T^, < T < T* with 7* = ep~^ 
and T*{i)=e-^'. 

Proof. As in subsection 2.2.1 consider shift with respect to ,^2 and it will 
be '^ iT. So, uncertainty principle requests £T x £ > h^^^ which is exactly 
our restriction to £ as T = ep"^ . Meanwhile xi < Cop^^ + Cp^^iT and it 
remains less than Cp~^ as 7~ < 1. On the other hand shift with respect to 
X2 does not exceed Cp^^ + Cp 2 T in virtue of proposition 2.17 below and it 
remains less than ei as T = ei^^^ unless p > i^ which is impossible in the 
transitional zone. 

We surely need to keep pi > h^~^ but one can check easily that pi > 

p-\h\-^ > h^-\ n 
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Corollary 2.16. (i) Contribution of the part of the transitional zone (1.14) 
where (2.38) is fulfilled to the remainder does not exceed C^^^h^^ ; 

(a) Contribution of the transitional zone (1.14) to the remainder does not 
exceed 

(2.39) Cfi-^h-'^+ 

C{fih)h-^-^mes9x{{x G dX, \VaxV\ < C(/i/7)5''^})/7-'~^ 
(Hi) In particular, it does not exceed 

(2.40) Cfi~H-^ + C{fih)h-^-\ 
provided 

{2Alf \VaxVF-'\ < e =^ ±Vl^VF-^ > e 

and expression (2.41) does not exceed Cfx^^h^^ as ft < /)^i3+'^. 

While either sign in (2.41) assures remainder estimate (2.40), in the 
future we will need to distinguish between different signs in this condition 
as dynamics will be different (see example 1.4). 

We finish this subsubsection by 

Proposition 2.17. An "average" propagation speed with respect to X2 in 

'^trans doCS not CXCCcd Cp^. 

Proof. Proof coincides with one of proposition 2.9 and corollary 2.10. 

Surely trajectories in the squeezed reflection zone are not transversal 
to boundary anymore after rescaling but we do not need it as instead of 
/i' 2 "'5- vicinity we can take e-vicinity and appeal f.e. to results of section 3.4 

[Ivr2]. 

2.2.4 Analysis in the inner zone 

Consider now inner zone deflned by (1.16). As we know, in this zone we 
must assume some non-degeneracy condition to get remainder estimate 
better than 0{fj,h~^). In Chapter 13 [Ivr2] we used few of such conditions 
rendering different remainder estimates; all of them boiled down to the 
remainder estimate 0{fi~^h~^) under different restrictions to /i. 
Consider first propagation of singularities. 
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Figure 6: Squeezed propagation with tangency, reflections or miss the 
boundary 



Proposition 2.18. Let Q G Sh,pi,p2,-yi.'y2 '^^'^ ^^^ ^ ^^ {li,l2, Pi, Pi) -vicinity 
(with respect to (x, ^)) of its symbol. Then long as trajectories started from 
Q remain in A'jnn defined by (1.16), (2.17) 

(2.42) U % = 6/° 'Qy 

where U^ is a propagator for operator A considered in X° the fixed vicinity 
of X (i.e. in the whole plane effectively). 

Proof. Proof is trivial as 

(2.43) P>P-= max(C/i-\ /iU^-^) 

as then we can apply Fourier integral operators as described in Chapter 13 
[Ivr2]. 

As p is defined by (2.17) one needs to apply the same arguments as in 
the proof of proposition 2.9. Albeit in that proof trajectories "comfortably 
met" dX and now they "comfortably miss" it. D 

Remark 2.19. This proposition however does not imply that results of sec- 
tion 13.3 [Ivr2] are automatically valid in our case as in that section we used 
partition elements which could be rather large in some directions (sometimes 
as large as x 1 which could be replaced by h^ as we apply microlocal stan- 
dard uncertainty principle rather than logarithmic uncertainty principle). 
Now boundary prevents us from doing this. 

Still, however, as /i < h^^ with sufficiently small exponent 5 > Q these 
arguments of section 13.3 [Ivr2] remain automatically valid and leaving easy 
details to the reader arriving to 
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Proposition 2.20. Under non- degeneracy condition (13.3.57)^ [Ivr2j with 
V = 1 i.e. 

(2.44L E \D''{VF-')\>e 

l<\a\<m 

contribution of the inner zone A'lnn to the remainder does not exceed Cfi^^h^^ 
as 11 < h~^ with sufficiently small exponent S = S{m) > 0. 

Combining witli results in inner and transitional zones where the same 
contributions were derived even without non-degeneracy condition we con- 
clude that 

Theorem 2.21. Let conditions (0.1)-(0.4), (2.2) and (2.12) and non- 
degeneracy condition (2.44)^ be fulfilled on supp^ where ip{x) G ^^(M}). 

Then the remainder does not exceed Cfj,~^h~^ as fi < h~^ with sufficiently 
small exponent 6 = 6{m) > 0. 

Proof. As in the Chapter 13 (see (13.3.49)^) [Ivr2] the contribution of A'inn 
to the remainder does not exceed 

(2.45) Cfi-H-^ I T-^ dx2d^2 + Cfih^^ mes2(Q2) 

with the same definition (13.3.50) [Ivr2] of Q2', here T = e/i2p2 and c/^2 = dp 
and integral is bounded as p < cp. D 

In almost all other cases we however need to take remark 2.19 into 
account. Still there is an exception: arguments linked to evolution with 
respect to ^2 work in a bit larger zone, namely {(x, ^) : ^2 < Q} and thus 
in zone {x, xi < Cifi^^} (which is a definition of Xbound)- 

2.2.5 No-critical point case 

Assume temporarily that condition (2.8) is fulfilled. Then, as we mentioned, 
theorem 2.6 and proposition 2.15 cover this zone. So we need to consider 
zone 

(2.46) {x : Cop"^ < xi < e}. 
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However under condition (2.8) we can take x 1 scale with respect to 
X2 and consider drift in the time direction in which ^2 increases; then we 
remain in the zone in question for time T* = e'/i. Meanwhile under this 
condition shift with respect to ^2 is observable as |t| > T^ = efi~^. This 
leads us to to conclusion that 

(2.47) Under assumption (2.8) contribution of this (2.46) zone to the Taube- 
rian remainder does not exceed C/x""*"/?""^ 

(i.e. the same as the contribution of zone {x : xi < Cofi^^} albeit factor 
H~-^ now comes as T"*"-*- rather than as the width of the zone). 
So we arrive to 

Theorem 2.22. Let conditions (0.1) -(0.4), (2.2) and (2.12) and non- 
degeneracy condition (2.8) be fulfilled on supp-?/' where ip{x) G ^q°°(M^). 

Then contribution of zone {x : dist(x, 9X) < ei} to the remainder does 
not exceed Cji^^h"^ as fi < h^^^ with arbitrarily small exponent 6 > 0. 

Non-degenerate critical point case. Consider now case when non- 
degeneracy condition (2.41) is fulfilled. Then obviously contribution of zone 
{xi < CoyU^-*-} to the remainder does not exceed 

(2.48) Cfi-H-^ + C{fih)'2-^h-^ 
where the second term comes from zone 

(2.49) {x : xi < Cofi~\ IW^^l < Q max(^-\ (^/?)^-'^) } 

and this is not as good as (2.40) because p <^ 1. 
Now in the zone 

(2.50) {x : Cofi'^ <xi < ei, | I/I4 J > Q max(p-\ (/i/?)^-'^) } 

we can use T* = efii with i = e\ l/l/xj] and we can even take T* = i^~^ ; then 
contribution of this zone to the remainder does not exceed C^^^h^^. 
Meanwhile contribution of zone 

(2.51) {x : Qfi-^ <xi < ei, 1 1/14 J < Q max(/i-\ (p/?)^-'^) } 
to the Tauberian remainder under this condition does not exceed 

(2.52) Cfi-H-^ + Cfi{fihY2-^h-\ 
So we arrive to 
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Theorem 2.23. Let conditions (0.1)-(0.4), (2.2) and (2.12) and non- 
degeneracy condition (2.41) be fulfilled on supp-?/' where ip{x) G ^q°°(]R^). 

Then contribution of zone {x : dist(x, 5X) < ei} to the remainder does 
not exceed (2.52); in particular, it does not exceed Cfi^^h^^ as ft < h^~"i . 

This is exactly the same estimate as was derived in Chapter 13 [Ivr2] 
without boundary under condition 

(2.53) |V\/F-^| + |VVF-^| >e. 

Now let us consider derivatives of VF^^ with respect to xi which leads 
to the drift in the direction X2. 
Our goal is to prove 

Theorem 2.24. (i) Let conditions (0.1) -(0.4), (2.2) and (2.12) and non- 
degeneracy condition (2.7) be fulfilled on supp-?/' where ^(x) G ^q°°(]R^). 
Then the remainder does not exceed 

(2.54) Cfi-^-^ + Cfi^h~^ 

as II < h^^2 ; in particular remainder does not exceed Cfi^^h^^ as fi < h^^3. 

(ii) Additionally assume that condition (2.41) is also fulfilled on supp?/'. 
Then the remainder does not exceed 

(2.55) Cfi-^-^ + Cfih'^-^ 

as II < h^^2; in particular remainder does not exceed Cfi^^h^^ as fi < h^^T . 
Proof, (i) First assume that 

(2.56) P>P-= max{Cfi-\ i^h'^-^) . 

Then along the whole trajectory xi > fi~^p > ji^^hi"^ and using n~^h- 
Fourier integral operators we can reduce operator to the canonical form 
of section 13.2 [Ivr2]. We can then notice that the shift with respect to 
"new" xi is fi^^Tk where k x IVVF^-*"! but the scale with respect to the dual 
variable is p~^p and we need to write uncertainty principle 



(2.57) fi-^kT X n-^p > fx-H 



li,l-5 
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or plugging T ^ ^~^ and /( x 1 we get 

(2.58) P>P-= max{Cfi-\ fi^h^'^) 

and this restriction is stronger than (2.56) which in turn is stronger than 
(2.17). 

Therefore 

(2.59) Under condition (1.8) contribution to the remainder of A'inn defined 
by (1.16), (2.58) does not exceed C^~^h~^ . 

Meanwhile with this choice of p contribution to the remainder of A'trans 
defined by (1.14) does not exceed Cph~^ which is exactly the second term 
in (2.54). 

Statement (i) is proven. 

(ii) To prove (ii) we need to note that in A'trans shift with respect to ^2 is 
observable as 7" x /i^-^ and | Wx^l x i satisfies (2.38). Then the contributions 
to the remainder of all zones save {Atrans, \^X2\ ^ ^} are 0{fj,~^h^^) and 
under condition (2.41) the contribution of this zone does not exceed Cpih~^ 
which is exactly the second term in (2.55). D 

Remark 2.25. (i) Conditions (1.8) and (2.41) are fulfilled generically at the 
boundary. 

(ii) Obviously in the setting of theorem 2.24 the inner zone provides the 
worst contribution to the remainder and we will need to improve it in frames 
of condition (2.24) using the strong magnetic field approach there. 

2.3 From Tauberian to magnetic Weyl 
formula 

Now our goal is to pass from Tauberian with T = epT^ to magnetic Weyl 
formula and estimate remainder Ri^^. We also consider extended Weyl 
formula and estimate remainder V(^ 



00 • 



Theorem 2.26. Letip G "^"^(X) he a fixed function with a compact support 
contained in the small vicinity of dX and let conditions (0.1)-(0.4), 

(2.60) F > eo 
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and 

(2.61) V<-eo Vx €6(0,1); 

be fulfilled there. Further, let condition (2.6) be fulfilled. 
Then 

(i) Under condition (2.8) both R^ defined by (2.14) and R^^ defined by 
(2.15) do not exceed Cfx^^h^^ as ji < h^^^; 

(a) Under condition (2.44)^ both R^ and R^"^ do not exceed Cfi^^h^^ as 

(Hi) Under condition (2.41) both R^ and R^'^ do not exceed (2.52) as 
/x < h^~^; in particular remainder does not exceed Cfi^^h^^ as ji < h^^^ ; 

(iv) Under condition (2.7) both R^ and R^"^ do not exceed (2.54) as /i < 
h^~^; in particular remainder does not exceed Cfi^^h^^ as fi < h^^3 ; 

(v) Under conditions (2.7) and (2.41) both R^ and R^^ do not exceed (2.55) 
as ft < h^~^; in particular remainder does not exceed Cfi^^h^^ as ji < h^^"' . 

Proof, (a) To go from R""" to R^ is easy: using condition (2.61) (i.e. | \/| x 1) 
we can apply the standard results of Chapters 4 and 7 [Ivr2] after rescaling 

X (-)• yUX. 

Going from R^ to R^^ is more subtle. Note first that we can drop all 
terms with m > in (2.14). Therefore only surviving terms are those with 
h^^[fihY" with integration over X and h~^(fihY" with integration over dX. 

(b) Consider first terms with integration over X. 

Note first that under condition (2.8) their sum is equal to 

(2.62) / Ar^^(x, ^/j, r)V^(x) c/x 

Jx 

modulo 0[h^^{fihyy, actually we go in the opposite direction: from expres- 
sion (2.62) to its decomposition into powers oi h = ^h. 

Consider expression (2.62) under condition (1.8); again we can replace 
d{T — V — {2j + l)Ffih) by derivative with respect to xi (of high order and 
with smooth coefficient) of smooth ^^-function; so integrating by parts we 
again arrive to the sum of the terms in question plus the similar integral 
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over the boundary; however the latter contains at least one extra integration 
with respect to r so we get 



(2.63) 



{^Mhf"h-'J{T -V- (2j + l)F^/));0„,,(x) ds. 



where r > 1. If we replace in the latter term summation with respect to j 
by integration the error will not exceed C{fih)''~^^h~^ i.e. Cfi^ which is less 
than (2.54). 

On the other hand, under condition (2.41) the error in question will not 
exceed C{^h)''^^^h~^ which is less than (2.55). So, in (iiv) and (v) we also 
can replace terms with integration over dx into A/''^^, may be changing k'^q. 

The same arguments albeit without integration with respect to xi work 
under condition (2.41) alone; however we gain only factor {fih)2. 

The similar arguments work in frames of (ii) as well. 

(c) Now case (i) becomes the most complicated as many terms should be 
taken into account. We apply the following trick: consider the same operator 
albeit we replace D2 everywhere by D3: 

(2.64) A:=J2~Pjg'\xi:X2)h+V{x,,X2), 

j.k 

Pi = hDi, P2 = hD^-fiV2{x^,X2) 

Then it will affect in R^ only terms with m > 1 we do not care about. 
But the problem remains microhyperbolic in the variable X2 and therefore 

everything works as it should. We need to consider then shifts with respect 
to .^2 only and therefore only averaging with respect to X2 is needed. 
Note then that what we get instead of J e(x, x, t)iIj dx is 



(2.65) 



(2.M-/e(x„x,x.,6,r)^.dx.dx,d& 



where e(xi,yi;x2, ^2, t) is a Schwartz kernel for spectral projector for 1- 
dimensional operator. Here integration with respect to X3 is not needed. 
Neither is needed integration with respect to ,^2 as we pass from Tauberian 
expression with T = e/i~^ for e(...) to e(...) itself. Finally we change in 

(2.65) 6 by 6. 

Also we can replace everywhere (save V2) xi by while V2 we replace by 

(2.66) V2= V{0,X2) + {d,,V2){0,X2)x,; 
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it will not affect essential boundary terms. But then boundary terms in 
R^ together must match to a boundary term in R^"^ . This proves (i) 
completely. 

(d) Exactly the same arguments work for (ii)-(v); we do not use condition 
(2.7) at all and in A'bound defined in terms of (^3,X2) we consider shifts after 
hops with respect to X3, so we do not need to integrate over X3. Without 
condition (2.21) we use the trivial 0{ph^^) estimate for contribution of 

'^trans- 

Under condition (2.21) in A'trans we consider shifts with respect ^2 and 
again integration over X3 and .^2 is not needed. D 



3 Strong magnetic field 

In this section we consider a case of the strong magnetic field when the 
results of the previous section are not as sharp as we want (so the remainder 
is not 0{fi~^h~^)). As under different assumption it happens under different 
restrictions to /x, we consider separately different cases. 



3.1 Most non-degenerate case 

If condition (2.8) is fulfilled, we need to consider only the case of very strong 
magnetic field 

(3.1) h^~^<fi<h-^ 

and then operator is X2-microhyperbolic. Therefore as ^ = ip'ip'' with 
ijj' = ■^'(x2) G ^0°^ ^'^d ■?/'" = '?/'"(xi7~-'-) ^) with ■?/'" supported in |xi| < C07, 

(3.2) \Ft^h-WXT{t)r^Pu\ < Cp~^h~^T{jy, 

and then 

(3.3) \F,^,-i,XT{t)Vi,u\ < Cfi-'h-' 



^) Averaging with respect to xi is not needed at all. 
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as h < T < T* = eo where fi~^ comes as a measure of Xbound and therefore 

(3.4) R"^ < C^-H-\ 

Further advancing method of successive approximations with unperturbed 
operator^) 

(3.5) A:=Y,Pjg^'{y)Pj+V{y). 



J.k 



Pj = hDj - ^Vj{y) - f, J2{xj - yj){d, Vj){y) 



we see that the first term results in expression (2.65) of the magnitude 
Cfi~-^h^^ while any next term acquires factor h in the corresponding power 
and thus does not exceed the remainder estimate. 

So, under conditions (2.8) and (3.1) and indicated ip the remainder does 
not exceed Cn~^h~^ while the principal part is given by the Tauberian 
expression for the first term in the successive approximation method. 

On the other hand, if we take ip" G ^o^(^^^) (supported in (|, 1)) and 
7 > /i^-*- we can take 

(3.6) T* X min ((rS)'^", f^l'^'") ■ 
Really, we could take 7"* x 7^"°" as 

(3.7) 7>7:= Co/i-' + Co/i-U^'^ 

which completely covers the case /i < h^^^. Otherwise we can arrive to this 
casescahngx i-)- x^, /x i-)- /x^, h i-)- h(~^ with ^ = min(l, {fj,^^h)2, [fi~'^ Iry'^Y) 
with large s. 

Therefore the contribution of the strip {x, xi x 7} to the remainder 
does not exceed 

C7r^((/)-S)-^-'^ + /i-S-^+-) 

and hence contribution of {xi < e} does not exceed this expression integrated 
over 7~-^a'7 resulting in 



C/,-^((/,-S)-^-'^7|,=^-: +/i-S'^l7=l) ^ Cf,-'h 



lu-1 



''' In comparison with (2.65), (2.66) we freeze at y not (0,y2) at this moment. 
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So, we can take ip fixed function ratlier tlian scaled witli respect to xi. 

However let us partition it into functions supported in {xi < 27} and 
in {xi > 7}. Then such expression in the latter case is not affected by 
the presence of the boundary resulting in the same expression but with 
approximation term calculated for the whole space. 

However in the former case the presence of the boundary should be taken 
into account. Let us use again the method of successive approximation but 
use as unperturbed operator one with xi set to everywhere save in the 
linear part of magnetic field; thus unperturbed operator is 

(3.8) A := h^Dl + (/ixi - hDjf - 1/1/(0, X2). 

Again as the main part of asymptotics is of magnitude 7/?"^ and each next 
term acquires factor 7, so only first two terms need to be considered. 

So, let us consider the second term; we claim that calculating this term 
one does not need to take a boundary into account. Really, as perturbation 
vanishes at the boundary one needs to kill xi before restricting to the 
boundary, but it can be done only by commutator and then factor h rather 
than 7 appears. It is not enough but if we plug instead of tp'J^ function ^" 
with Cofi~^ < 7 < 7 and supp^^" disjoint from we acquire factor {fi'~f)~^ 
and then contribution to the error is C'~fh~^{fi'~f)~^ and it boils down to 
C/i~"^ /?""'" after summation with respect to 7. 

Then we get the final answer as the sum of two terms: one is for operator 
(3.5) albeit with calculation (before taking f) in the whole plane i.e. 



(3.9) 



/)-2 f AT^W^^^ /i/?)^(x) dx 

Jx 



and the second one for operator (3.8) but in half-plane and subtracting the 
same expression for the same operator albeit in the whole plane we arrive to 

(3.10) J (y ei(xi, xi; X2, 6, 0) c/^ - h-^^f^'^{x, /i/)))^(x) dx 

where ei(xi, yi; X2, .^2. 0) is the Schwartz kernel of one-dimensional operator 

(3.11) /J^D^ + (/iXi- 6)2- 1/1/(0, X2) 

and Aq^^ refers to TV''^^ calculated for the same operator (3.8). 

We will transform operators (3.8), (3.11) and silently (without changing 
notations) transform e(...) and A/g^^. 
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Changing Xi ^-^ Xi( and ,^2 ^-^ ^.iCf^ we acquire factor fi( and get instead 
of the first term in (3.10) 

(3.12) {2nh)-\fi j e,{x,.x,:x2.^2.0)^'{x2)r{xiOdx,dx2d^2 
meanwhile transforming operator into 

(3.13) h'c' [of + /)- VC'(xi - 6)' - W{0, X2)h-'c 



-2a-2 



One can drop a factor in the front of operator and select ( = fi 2/12 thus 
resuhing in the answer 

(3.14) {2n)-^fih-'2 / ei(xi,xi;x2,6.0)^'(x2)V'"(xi/i5/,-5)c/xic/x2c/6 

and in operator 

(3.15) D2 + (xi - 6)' - l^'^h-^ 1/1/(0, X2). 

However we need to subtract from (3.14) also transformed the second term 
in (3.10). We can then tend ( — ?■ +00 (the error will be negligible) and the 
total difference will tend to 

(3.16) h' [ ^Tou..{x2.^^h)i^{x)ds,. 

JdX 

So, the final answer is 

(3.17) h-' f Af''^{x,f^h)^{x)dx+h-' I ACZnd(^2,^/')^(x)c/s, 

Jx JdX 

with * = D, N and we arrive to 

Theorem 3.1. Let ip G ^'^{X) be a fixed function with a compact support 
contained in the small vicinity of dX and let conditions (0.1)-(0.4), (2.60) 
and non- degeneracy condition (2.8) be fulfilled on supp-?/'. 

Then formula (3.17) gives N with an error R'^^ = 0{fj,^^h^^) as 11 < h^^ . 
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3.2 Generic case. Analysis in inner zone 

Now we are interested to improve results of the previous section in the 
generic case i.e. when both conditions (2.7) and (2.41) are fulfilled. Then in 
virtue of theorem 2.26(v) we can assume that 



(3.18) fi>h^'y. 

We start from the simpler analysis in A'lnn- In this case as p > p = 
{fj,h)3h~^ we need to consider operator without boundary condition; however 
presence of the boundary as we remember manifests itself through uncertainty 
principle; we should take 7"* = h^^^ i^'^ or 7"* = jih^^^ p^^ whatever is smaller 
and at this moment we are interested only in the zone where 7"* > ep~^. 

Actually we can take here effectively even 5' = in the following sense: 
note that 

(3.19) \F,^,-irMiY{^%)\ < Cfi^H-^pix (fiT + l) 
as0G<^o°^([-l,l]), r> /7and 

(3.20) \F,^,-.,XT{t)\-{U%)\ < Cp~H-'pix {pT + l){^y' 

as X e '^o°^([-l, -\] n [|, 1]) and T^ < T < T* with unspecified at this 
moment 7"* and 



(3.21) 7, = min(/i/)p-\r^)/? = /?x ^ ^_2 ^ ^ 



pp ^ as p > pi'^, 

i~^ as p < pf 



Recall that Q = Q{x2, hD2) is an (£, p) admissible element with pi > h^~^. 

Really for one winding (i.e. T = ep~^) estimate (3.19) is obvious, and 
we need to take in account A/ x /i7" + 1 windings. 

Estimates (3.19) and (3.20) imply that 

(3.22) \Ft^i,-ir(j)T{t)V{U%) I < Cp-^h-^p£ X {pT, + l) 

as G '^o°^([-l, 1]) and 7", < 7 < T*. 

Remark 3.2. Therefore from the point of view of the remainder estimate, 
rather than the final formula we need to take in account A/* x (/iT,, + 1) 
windings but in the main part of asymptotics still {pT^:h~^ + 1) windings 
should be taken in account. 
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Let us discuss T*. As average propagation speed with respect to x does 
not exceed Cqji^^ and propagation speed with respect to ^2 does not exceed 
Co^, our dynamics remains in the same (£, p)-element for time 

(3.23) 7* = emin(/i£,/)r^) 

and 7; : 7* X hp-H~\ 

Therefore contribution of [i, p)-element to the Tauberian remainder with 
T > {T^h^^ + efi^^) does not exceed 

(3.24) Cfi-^h-^p£{i2T, + l) X T*-\ 

Note that in zone i < /i^2p2 we can reset i = i{p) := pT^pi thus covering 
the whole zone with a fixed magnitude of p by a single element. Also note 
that (!,{p)p = p~2p2 > p2h > h^~^ as p > p. 
Therefore 

(3.25) Contribution of (£(p) , p)-element to the Tauberian remainder with 
T = ph^~^ p~^ does not exceed 

Cp-'h-'i{pY{^ + 1) X c + cp-^-'p 

and summation with respect to p G [p, ep] results in 

(3.26) C\\ogh\ + Cp-H-\ 
So, 

(3.27) Contribution of zone {p > p, pi^ < p} to the Tauberian remainder 
with T = ph^~^ p~^ does not exceed (3.26). 

Remark 3.3. We need to keep i{p) > Cp~^ but this is definitely case as 

2 

p > h^5 and we are ensured in this by condition (3.18). 

Similarly, consider case pi'^ > p > P- Then automatically i > i{p) and 
ip > h^'^ and 

(3.28) In this zone we can reset p = p(£) = p(?. 
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Really, to avoid collision with A'trans "we select time direction in which p (and 
thus .^2) increases. It is possible because as long as |t| < efii we remain in 
the same e£- vicinity. Again i > Cfi~^ in virtue of condition (3.18). 
Then 



(3.29) Contribution of {i, p(£))-element to the Tauberian remainder with 
T = h^~^ i~^ does not exceed 



./i/? 



C^i~'h-'p{i){^ + l)>,C+Cp~'h 



lu-lfl2 



and summation with respect to £ G [£, e] results in (3.26) where i = fi 2p2. 
So, 

(3.30) Contribution of zone {p > p, pi^ > p} to the Tauberian remainder 
with T = ph^~^ p~^ does not exceed (3.26). 



P fd' < p 


/ 








^if>P 




"Abound U A'trans 



fp 




(a) fif < p (b) ^if > p 

Figure 7: Partition of two zones in A'inn 



pf > P 



So we arrive to 

Proposition 3.4. Under assumptions (2.41) and (2.7) contribution of X\„n 
to the Tauberian remainder with T = ph^~^ p^^ does not exceed (3.26). 

We are completely happy with this estimate unless p > h^^\ log h\^^ and 
we need to derive contribution of much more troublesome zone A'trans before 
thinking if we should improve it. 
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3.3 Generic case. Analysis in transitional 
zone 

Let us consider zone A'trans- Recall that its contribution to the Tauberian 
remainder with T = ejj,~^ (and thus with T = h^^^) does not exceed 
Cph~^ = C{fih)3h^^^^ in the general case and Cpih~^ = C{fih)6h^^~^ under 
condition (2.41) and therefore (as we assume (2.41)) we should consider only 
case 

(3.31) fx>h^-T3. 

Definitely A'trans is leaner than X\„„ (albeit as /i > h^~^ it is thick enough 
to eliminate A'bound)- However the main problem there is that in the current 
settings we are not aware of any lower bound of the propagation speed and 
the only bound we know is an upper bound Qpz in both directions. 

Therefore there is no mechanism except drift with respect to ,^2 to break 
periodicity and we must take 

(3.32) 7, = ehe-^ 
and 

(3.33) T* = ep-k. 



Sure we need to have T^, < T*, i.e. i > io with 
(3.34) £0 :=p5/,5 =(^/,)i/,5-^ 

Obviously io < i = {ph)ih~^ unless p < h^7 the case we are not interested 



m. 



Remark 3.5. Note that it it requires time t* x pi~^ to pierce through A'trans 
and t* < T* a.s i < i which is the case we are interested in. Otherwise we 
would be able to increase 7"* further. 

So, in the same manner as before contribution of {i, p)-element with 
i > Iq to the Tauberian remainder with any T > T^h~^ ^^ by 

(3.35) Cp-^pih-\pT, + 1) X T*-^ = Cph-^ + Cp-H-^p'2 

8) As 7* /?-■*' > T* we reset T ^ T* . 
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and summation over i > io results in 

(3.36) Q542 ^ C/i"i/,-ipf I |og/,| X C(/i/?) 5 /j-i-^ 

On the other hand, contribution of {£0 , p)-e\em.ent to the Tauberian remain- 
der with T > eji~^ does not exceed 

(3.37) Cploh-^ 

which coincides with the (2.37). 
So we arrive to 

Proposition 3.6. Under condition (2.41) contribution o/Atrans with any 
T > T^h^^ = h^~H^^ (reset to T > h^'^ as either £ < Iq or £ > I; see also 
^^) does not exceed (3.36). 

In particular, it does not exceed Cfi^^h^^ as ji < h^^s. 

Remark 3.7. (i) Surely we need to keep £ > Cqji^^ thus requiring £q > Cqjj,^^ 

2 
i.e. jj, > h^5 but this is a case. 

(ii) We also need to keep an upper bound to speed greater than Co/i^"*", i.e. 
p2 > Co/i~"^ i.e. fi > h~'i but this is also the case. 

(iii) Further, we need to keep £op > h^~^ i.e. p> h^^ , but this is again the 
case. 



'^trans 

Figure 8: Zone A'trans 

Can we increase T* in these arguments? We need to do it only under 
condition (3.18). Then as we will prove later 

Proposition 3.8. Let 

(3.38) Dirichlet boundary condition be given on dX . 
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Then in the transitional zone {\p\ < p} as 

(3.39) p > h--^-^ 

average propagation speed in one direction (direction of hops when we have 
chosen time) is bounded by Cp^ and in the opposite direction it is bounded 
by Cp-^. 

Note, that condition (3.39) ensures that ps ^ p~'^. 

So, assume that (3.38) holds. Then selecting the time direction for given 
partition element so that i increases along hops we can chose 

(3.40) T* X 1. 

Really, we need to leave ^-element and come back. This comeback includes 
either going in the direction opposite of hops with the speed not exceeding 
Cqp~^ and requires time x pi which is larger than (3.40) as £ > p~-^ and 
/i, > /?" 4 or leaving zone i < e which requires at least time x 1 ^) . 

Actually T* x 1 is better (larger) than given by (3.33) we used before 
only as £ < p2, but this is only zone we need to care about. 

So, contribution of the given element to the remainder does not exceed 
the left-hand expression of (3.35) which now becomes 

(3.41) Cp-^pih-^ X {phr^ + 1) X m\n{pk-\ l) x Cpr^ + Cp-^h-^pi. 

Summation of the second term in the right-hand expression over partition 
results in 0{p~^h~^). Summation of the first term in the right-hand ex- 
pression over partition with i > ii results in the same term calculated as 
i = ii and coincides with Cpiih^^ which also estimates the contribution of 
zone {£ < i^}. Here again we chose ij from condition 7",, < T* i.e. now we 
replace Eq defined by (3.34) by 

(3.42) h := h'2 < io 
resulting in 

(3.43) Cph-'2 = Cph's^^. 

Therefore (with pending proposition 3.8) we arrive 

^' As X is a bounded domain and F x 1 in X this would mean run-around X along 
boundary which as we see later is impossible under certain assumptions. 
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Proposition 3.9. Assume that Dirichlet boundary condition is given on 
dX . Then under assumption (2.40) the total contribution o/ A'trans to the 
Tauberian remainder does not exceed (3.43). 

In particular, it does not exceed Cfi^^h^^ as fj, < h^^w and it is always 
has an extra factor fii h6~^ in comparison with h~^. 

Can we do better than this? Yes, as we can take 

(3.44) T* X /i£. 

provided ^2 and i increase in the same time direction (on given ^-element) 
as F12 < 0. 

Really, then selecting such time direction we remain in A'trans H X\r,n and 
as we reach £ x 1 dynamics is already in X\„r, where speed in both directions 
is Oifi-^). 

Assuming that the critical point of W := VF^-^\dx is X2 = we note that 
hops are to the left (as F12 < 0). So, we select time direction of the sign of 
(— X2), so £ X |x2| increases. 

Meanwhile those expressions have the same sign: d^2/dt, —d^^ 1/1/ and 
(—9^2 l/l/)x2 and our condition means that these expressions are positive, 
i.e. condition (2.41)~ is fulfilled which matches to cases displayed at Fig- 
ure 4(b), (d). So, under this condition we can select T* according to (3.44), 
therefore replacing li defined by (3.42) by 



1.1 ^ 



(3.45) £2 :=/U"5/,3-'^ < 



'1 



(as fJ, > h 2 '^) from condition 7"* < 7"* and arrive to remainder estimate 
(3.46) Cphh-^ + Cfi-^-^ = Cfih-^ + Cfi'H~\ 

So we have proven 

Proposition 3.10. Assume that Dirichlet boundary condition is given on 
dX . Then under assumption (2.41)~ the total contribution of Xt^ans to the 
Tauberian remainder does not exceed (3.46). 

In particular, it does not exceed Cji^^h"^ as fi < h^^^ and it has an 
extra factor fi^ h^^^ in comparison with h^^ (as h^~^ < /U < h~^). 
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Finally, assume that hops and magnetic drift have the same direction 
near point in question which is the case under condition 

(2.7)+ \d.,VF~'\ < e ^ d,,VF~' > eol 

compare with original (2.7). 

Then according to proposition 3.8 there can be no roll-back, only run- 
around with the path of the length x 1. This corresponds to Figure 4(a), (b). 

However if condition (2.41)^ is also fulfilled (so we are in frames of 
Figure 4(b)) fast run-around along boundary is not possible either, so the 
run around must contain the segment of the length x 1 inside X\nn (even 
inside {x, xi > e}) but the speed here is 0{fi~^) and therefore T* x fi. Then 
L, < T* defines 



(3.47) i^ ■= ^—2 h-2-' <^ 



'-2, 



this leads to remainder estimate Cpish ^ = n^h^ ^ which is 0(;U ^h ^) as 
/i < h^~^ and 0{h~^) otherwise. Thus we have proven 

Proposition 3.11. Assume that Dirichlet boundary condition is given on 
dX. Then under assumptions (2.41)^ and (2.7) the total contribution of 
'^trans to the Tauberiau remainder does not exceed Cfi^^h"^ as ft < h^^^ and 
Ch^^ otherwise. 

Conjecture 3.12. The above Tauberian estimates hold with 5 = 0. 



To prove this conjecture we must reload A'bound and X\nn and to derive 
there the same remainder estimates as we proved in X\„„ but with p = {ph)3. 
It does not look difficult as we need just to rescale x, — )■ X7 with 7 = h[p/p)~^ 
with arbitrarily large k. However then h \-^ h = {p/p)~^ and it is not small 
enough to keep number of reflections below h~^ . Really, we need to consider 
at least phi~^ number of reflections which means phi~^ < h^^ and as h goes 
to 1 we need to allow i > {ph)2 only! 

Another approach would be based on logarithmic uncertainty principle 
(albeit instead of h~^ logarithmic factor | \ogh\' would appear) but we do 
not have here theory similar to one of section 2.3 [Ivr2]. 



Chapter 3. Strong magnetic Geld 46 

3.4 Propagation of singularities in 
transitional zone 

Let us study propagation of singularities in A'trans; later we assume that the 
Dirichlet boundary condition is given. We will use the technique developed 
section 3.4 [Ivr2]. However there is a large difference: in section 3.4 [Ivr2] 
we treated basically trajectories with a single tangent point while here we 
need to treat fiT such points in one shot. 

Consider real function z/''(x, ^, t, r) and by Weierstrass theorem replace 
it modulo (r — a(x, ^)) by a linear with respect to ^i function 

(3.48) ^'{x, e, t, r) = Mx. 6, t, t) + ^i(x, 6, t, r)^ = 

ipix, C t, r) + a{x, ^, t, r) (r - a{x, 0) ■ 

Then as coefficient at C,i in a(x, ^) is 1 

(3.49) {r -a,tlj} = {T- a, ^} + ai (r - a(x, 0) 

is at most quadratic with respect to C,i and therefore the following identities 
hold: 

{{hDt - A)u, ^™ty) = {u, {hDt - A)tlj''u) + ih{hDiu, tp'"u)9x- 

ih{u,hD^ru)dx. 
{^""u.ihDt - A)u) = {u.ij'^ihDt - A)u) + ihiij'^u.ihDt - A)u)ax 

implying identity 

(3.50) 2h-^ Re i{{hDt - A)u, ^™ty) = -{ih-^[{hDt - A), ^™]ty, u)- 

{hD^u, ij'"u)ax + {u, hD^tlj''u)9x + {tp^u, [hD, - A)u)ax. 

In this identity terms containing {hD^ — A)u will be either negligible or under 
our control anyway; {{hD^ — A), ip'*'] = —lh{T — a, ip}^ modulo smaller terms 
and only boundary terms 

(3.51) - {hD,u, V^"u)9x + {u, hD^^''u)9x 

need special analysis. Under Dirichlet boundary condition we rewrite them 
as 



(3.52) -{hDiu.ip'^hDiu) 



dX 
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and under Neumann boundary condition we rewrite them as 

(3.53) {u, h'Dlr^u)Qx = -{u, {hD, - /l)^r")ax + {u, {hD, - A')^P'^u) 

where A' = A — h^D^ and the difference is that in the former case definiteness 
of this quadratic form is ensured by ipi having definite sign but in the latter 
case not. 

Remark 3.13. (i) The same difference would manifest itself as we would 
try to prove results of section 3.4 [Ivr2] for wave equation under Dirichlet 
and Neumann boundary condition. The truth is that the former satisfies 
Lopatinski condition while the latter does not: uniformity breaks in the 
tangent zone. 

(ii) The same difference manifests itself through different behavior of eigen- 
values Adj(?7) and A|\ij(^) as ?7 — )■ +oo: while they both tend to {2j + l)/i/7 
the former are tending to it from above and the latter from below. This 
implies that under Neumann boundary condition at least for model operator 
some singularities propagate in the direction opposite to hops and this 
propagation is faster than magnetic drift. Such difference will be the most 
transparent in the case of very strong and superstrong magnetic field in 
section 4. 

So, we assume that Dirichlet boundary condition is given on dX. We 
select 

(3.54) ^1 > 
and thus we should request 

(3.55) {r - a, V^'} + /3(r - a) < with /3 = /3(x, t, ^2, r). 
We take originally 

(3.56) i'' = x{(t>), <P = <l^{x.O 

where '^°°(]R) 9 x is a function of the same type as in section 3.4 [Ivr2]: 
namely supported in (—00, 0] and with x' < on (— oo,0). 

Then as r — a(x, ^) = has two real roots ±ri 7^ we conclude that 

(3.57) ^1 = ^ (xmi = v)) - xmi = -V))) > 
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provided 

(3.58) - %0 > eo- 

It does not satisfy (3.54) as r — a' < but we will handle this in the same 
way as in section 3.4 [Ivr2]). 

Now, modulo terms of the type (3 ■ {r — a) 

(3.59) {r - a, ^} = {r - a, ^'} = -x'(0) ■ {r - a, 0} = 

-xl{4>){r-a,<P}^-A'{r-a,^} 

cts Xi = V~X' is a smooth function (under correct choice of x) and ip is 
again linear with respect to ^i function, 

(3.60) ^ = xi(0)({r-a,0})5 
and we used Weierstrass theorem again. We want 

(3.61) {r-a,0}>eo. 

Remark 3.14- Conditions (3.61) and (3.58) mean respectively that in- 
creases along trajectories as xi > and at reflections. 

Let us recall that instead of (3.54) we actually have a weaker inequality 
Namely instead of (3.54) we have 

(3.54)' iJ^ > -X2(0)'x3(r - a'f 

with both X2 and Xa smooth functions supported in (— oo, 0] and notice that 
both operator {hDt — A) and Dirichlet boundary problem for it are elliptic 
as r — a' < and we can apply elliptic arguments there. 

So, repeating arguments of section 3.4 [Ivr2]) we conclude that if 

(3.62) \NF'+\{hDt - A)u) n {0+ < e} n {0 < t < 7} = 0, 

(3.63) WF'(ty|ax) n {0+ < e} n {0 < t < 7} = 0, 

(3.64) (WF'(ty) U \NF'{hDiu)) n {0+ < e} n {t = 0} = 
and 

(3.65) (WF'-"(ty) U \NF'-^hDiu)) n {0+ < e} n {0 < t < 7} = 0, 

(3.66) WF'-"(/?Dity|ax) n {0+ < &} n {0 < t < 7} = 
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then 

(3.67) (WF^(u) U \NF'{hDiu)) n {0+ < 0} n {0 < t < 7} = 0, 

(3.68) \Nf'{hDiu\9x) n {0+ < 0} n {0 < t < 7} = 0. 

Here a > is a sufficiently small exponents, e > is an arbitrarily small 
constant and 7 > is an arbitrary constant, 

(3.69) 0^(xi,X2,6,t,T) = 0|5,^±„ r] = {T-a'y2 

and WF^ are defined in terms of pseudo-differential operators b{x, t, hD2, hDt). 

As we can plug (0 — e) instead of we by induction can get rid off 
assumptions (3.65), (3.66) (assuming that u is temperate). 

We also can rescale x i— )■ 7x (and h i— )■ hT~^, fi \-^ fiT) thus replacing e 
by e7 when we come back and we can also consider large 7. 

As a = ^1 + (.^2 — fJ'Xi) + V we select 

(3.70) = /tt + (/iX2-ei). 
In more general case we select 

(3.71) = /tt + /iX2 + F12P1. 
Note that (3.49) is fulfilled and also 

(3.72) {T-a,(t)} = K-d^,V; 

and (3.61) becomes k, — dx^V > e; so is fulfilled in two cases: 

(3.73) K = Co 
and 

(3.74) dx,V>eo. |fi:| < ei 

which allows us to prove respectively that "(//X2 — .^1)" propagates with a 
speed bounded from below by — Q and, as d^^ V' > eo by ei. To get rid off 
^1 we must to pass from 0+ to but we need to notice that 0"'' < + Cq. 
Then we arrive to 

Theorem 3.15. Let Dirichlet boundary condition he given on dX and let 
F12 < i.e. hops go to the left. Let x e '^o°^([|, 1]). Then 
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(i) For Co<T<T* 

(3.75) \Ft^h-wXT{t)Mx2)UMy2)\ < Ch^ 
as 

(3.76) X2 < —Cojj,^^T + y2 Vx2 G supp-?/'2, Vy2 G supp-?/'i; 

(^MJ Under condition (3.47) for Cq < T < T* (3.75) holds as 

(3.77) ^2 < eoAi "^ "?" + y2 Vx2 G supp-(/'2, Vy2 G supp-?/'i. 

This theorem imphes proposition 3.8 and thus justifies the improved 
results of the previous subsection. 

3.5 Calculations 

Now we need to move from Tauberian to more exphcit expressions for the 
principal part of asymptotics. 

First of all recall that the Tauberian expression for r(e(., ., 0) ^Qy) is 

(3.78) h'f F,^,-^,XTit)r{UVy)dT 

J — oo 

and if we replace Xr^t) by XT[t) ^°^ we will get 

(3.79) T-^ j F,^,-,,XT{t)r{U%)dT 

where x{t) = it~^x{t)- 

Let us list different cases depending on our assumptions 

3.5.1 General discussion 

As Q is supported in A'trans an absolute value of this expression does not 
exceed 

(3.80) Cfi-^peh-^ T-^ X fiT= Cph-^ 

10) Where as usual x> X "= '^^o°{[^^' 1]) ^-nd x = 1, x = on [— j, j] and later we will 
decompose xr in the sum of Xfi^) ^'^'^ X2-"r(*') 
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as 7~ > cq/U^-^ and therefore an approximation error should not exceed 

(3.81) Cpih-^ X Th-^A 

where A is the size of perturbation and we know that as we replace \/(xi, X2) 
by V(xi,y2) 

(3.82) A = Coi{8x2) + Co(6x2)2 = Qiifi-' + v7) + Co(^-^ + yTf 

as 6x2 = (/i^"*" + v7") and in the first term factor i comes as bound for \dx2 V\. 
Here v x p2 is an upper bound for propagation speed in A'trans- However we 
will justify that under Dirichlet boundary condition we can take in these 
calculations v x fj,~^. The difference between ps = (p/))3~'^ increases as p 
increases. 

Neglecting the second term in (3.82) (to be justified later), we acquire 
factor 

(3.83) Th-^A = a{fi-^ + V 7) Th-^ 
which after plugging T = h^^'^i^^ becomes 

(3.84) Cifi'^ Th-^ + Cv 72/j-i = C U~H~^ + Cyhr^) h'" 
and it is 0{h") as 

(3.85) £>£:= max(/i-\v5/?5)/)"'' 
and the contribution of zone l < d. does not exceed 

(3.86) Cpih^^ = Cpmax{fi-\yh'3)h-^-'' 

and as p < h^^^ modulo 0{fi^^h^^) it is 

a^lT^ as V X p2 
p^h " as V X p "^ 

which is 0{fi~^h~^) as either v x p2 and p < /)^~i3 or v x p^^ and p < /7'^~4. 
Meanwhile contribution of an element with i > £ does not exceed (3.81) 
which is 

(3.88) Cpih-^ X U~H-^ + CvhtAh-" 
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and summation over £ > £ results in the same expression with £ = £ and it 
is the same (3.87) as before. 

On the other hand, (the ratio of) second term in (3.82) to the first one 
is (/i^"*^ + vT) to £ and £ > ^~^ so only comparison is yh£~^ vs £ and £ is 
larger by the choice. 

So, 

(3.89) Under condition (2.41) contribution of A'trans to approximation error 
is 0(/X9/7 9^'^) for Neumann boundary condition and 0(^nih^^^ for Dirichlet 
boundary condition. 

3.5.2 Case of condition (2.7)^ 

Assume now that Dirichlet boundary condition is given and condition (2.7) 
is fulfilled. Then 

(3.90) r. = /,-'^min(M,^) 

and therefore the above arguments should be applied only as £ > £i with 

(3.91) £i=/i-5/,5; 

otherwise we replace £ by £\ (single element) and T = fi^hi^"^ and 5x2 = 
6xi = 0{fj,~^); so A X £ifj,~^. Here we assume that £i > Cjj,~^ as otherwise 
we reset it to this value which would lead to the error less than Cp~^h~^~'^ . 
Then obviously the total error under the same replacement V(xi, X2) by 
^(0,3/2) + (<9xi V')(0,y2)xi does not exceed 

(3.92) Cphh-^ X {fih)h-^-''hfi-^ X d^h-^-'^ = Z)-'^ min(^, ^) 

P £1 

which in turn is 0(^(i^^h^^ + h^^Y 

3.5.3 Non-degenerate case 

Finally, under condition (2.7) we take T^, = h^^'^ and A = pT^h^'' so each 
next term acquires factor Cp^^h^" i.e. the second term is Cph^^ x pT^h^'^ 
which is 0{ji^^ h^^) unless ji > h^^^, in which case it is 0{h^^). However in 
the latter case the the third term is 0(1) so we will need to consider the 
second term in approximations as well. 
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3.5.4 Boundary zone 

Here we take 7"^, = h^^"^ and A = fi'^h"'^ so each next term acquires factor 
C^^^h^" i.e. the second term is Ch^^ x ^^^h^" which is 0{^~^h^^^"). 
However the third term is 0(1) so we wiU need to consider the second term 
in approximations as weU. 

3.5.5 Inner zone 

In inner zone we make usual transformation and apply the standard method 
of successive approximations at point y rather than (0,72). 

3.5.6 Calculations. I 

As a result in the zone {xi > Cq^^^} we get 



h~^ fAf'''^{x,T,fih)^P{x)dx 



and in the zone {xi < Cofi'^} we get something like this but with J\f^^ 
replaced by A/'J^^ temporarily denoting eigenvalue counting function for 
operator with the boundary conditions. So this main part would come from 
parametrix G^ for operator hDf — Am the direction of ±t > (see previous 
and similar chapters) and this parametrix is equal to G^ + G^ where "bar" 
refers to freezing coefficients in y. 

Consider first term with Gq . While Gq is what it was on M^, V is not 
as now old V is replaced by r9(xi) = V — r0(— xi) and only the first term 
results in /)-W^^(y,r,/i/?). 

Note that contribution of terms with rQ{—xi)GQ and VQ{xi)G^ into 
main part of asymptotics are of magnitudes ix~^h~^ x nh = h~^ where /i""*" 
is the width of Xbound and fih is a semiclassical parameter after rescaling. 
Therefore replacing yi by will bring really the estimate we had referred to 
and these terms together will result in h'^N'^'^ 



corr 



3.5.7 Justification 

To justify setting of the upper bound of propagation speed fi~^ rather than 
p2 under Dirichlet boundary conditions we just note that if A' = A — A' then 
considering sandwich G^A'G^ ... A'G^S (one of G^ may be replaced by G^) 
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we can decompose A' = Aq + A'_^ + A'_ where A'q, A'_^_ and A'_ are copies of 
operator A' localized as |x2 —72! < 2co/i~"^(7" + 1), (x2 — 72) > |c'o/i^"^(7" + 1) 
and (x2 — 3/2) < — |co/i~"^(7" + 1) respectively. 

Note that due to the fact that in direction of ±t propagation speed does 
not exceed Cofi^^ we conclude that as we consider time direction ±t > 
sandwiches containing at least one factor A'^ are negligible, and sandwiches 
containing at least one factor A'-^ become negligible after we apply X2 = y2- 

So, the first term in approximation results in what we claimed as an 
answer. 

3.5.8 Calculations. II 

Meanwhile as we replace V(xi,y2) by V(0,y2) + (<9xi V)(0,y2)xi we notice 
that A = Co(6xi)^ x Co/i~^ as 5xi = 0{fi^^) which is smaller than (3.82). 

The same arguments show that removing (<9xi V)(0,y2)xi leads to the 
error not exceeding 

(3.93) C/i"5p/,-i = Cfxh-'^-"" 

as we need to use ii = fj,~2 instead of i defined by (3.86). 

Unfortunately expression (3.93) may be larger than anything we got 

1 

before. Fortunately there are alternatives: first, replacing xi by Wq (X2) 
leads to an error in operator 0{n~^p) and to the approximation error 

(3.94) Cifi-^pf^ph-^ = Cp^h-"; 

this would mean replacing U^ZnA^^) by ACTund(>^i) with W, = W^ + 

Second, we can replace xi by ^'^^2 i-e. use yV[iJi~^^2,X2) instead of 
lVo(-'^2) and we do not need to remove linear term [d^^ ^)(0. y2){x-i — f^'^^i) 
but just simply shift with respect to ,^2- This would lead to 0{fj,~^) shift in 
W[fi~^C,2,X2) and this shift could be ignored. 

3.5.9 Calculations. Ill 

Consider the second term in approximations. There are two cases when we 
need to do this: in the boundary zone as /i < h^~^ and under condition (2.8) 
as fi> h^-^. 
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In the former case however we need to consider only T < efi and then 
only one winding should be considered and perturbation does not exceed 
Co/i^"^, and we can take T = h^~^ so in fact we get an extra factor n~^h~". 
However scaling shows that in fact as we replace xt by Xt the error will be 
Cfi^^Th^^{h/Ty and summation results in Cfi^^h^^ error. Alternatively 
we could employ semiclassical approximation to show the same. 

In the latter case there is no boundary zone anymore. However scaling 
shows that the same arguments work (without going to semiclassics). 

3.6 Final results 

Let us compare Tauberian remainder and approximation error. 

(i) Under assumptions (2.8) both Tauberian remainder and an approximation 
error are 0{fi~^h~^). 

(ii) Under assumptions (2.7) and (2.41) Tauberian remainder is estimated by 
(3.36) + 0{fi^^h^^) and an approximation error is estimated by (3.87)+(3.93) 
and estimates (3.36) and (3.87) coincide, so the total error is (3.87) + (3.93) + 
0(/i-i/?-i). 

(iii) Under assumptions (2.7) and (2.41) and Dirichlet boundary condition 
Tauberian remainder is estimated by (3.43) + 0{fi~^h~^) and approximation 
error is estimated by (3.87) +(3.93) and (3.43) is larger than (3.87) so the 
total error is (3.43) + (3.93) + 0{i2-^h-^). 

(iv) Under assumptions(2.7) and (2.41)~ and Dirichlet boundary condition 
Tauberian remainder is estimated by (3.46) and approximation error is 
(3.93) + 0{fi~^h-^ + h~^) and the total error is (3.46)+(3.93). 

(v) Under assumptions (2.7) and (2.41)^ and Dirichlet boundary condition 
Tauberian remainder is 0{fj,~^h~^ + h~^) and an approximation error is 
(3.93) + 0{fi~^h~^ + h~^) and this is a total error. 
Thus we arrive to 

Theorem 3.16. Letip G ^°°(X) be a fixed function with a compact support 
contained in the small vicinity of dX and let conditions (2.60)-(0.4), (2.60) 
and and (2.61) be fulfilled on supp-;/'. 

Further, let condition (3.1) be fulfilled. Then 

(i) Under non- degeneracy condition (2.8) R'^^ is 0{fi~^h~^); 
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(ii) Under non- degeneracy conditions (2.7) and (2.41) Ri^^ does not exceed 
(3.87) + (3.93) + 0{i2-^h-^); 

(Hi) Under non- degeneracy conditions (2.7) and (2.41) R^^ does not exceed 
(3.43) + (3.93) + Oifi-^h-^); 

(iv) Under assumptions (2.7) and (2.41)~ R^^ does not exceed (3.46) + 
(3.93); 

(v) Under non- degeneracy conditions (2.7)^ and (2.41)~ Rq^ does not exceed 
(3.93) + 0(/i-i/j-i + /7-^). 

We leave to the reader 

Problem 3.17. Write down correction terms AAJ^b^un^j <-orr corresponding to 
the procedures described subsubsection "Calculations. II" . Then R^)!^ will 
be the same as in theorem 3.16 albeit without (3.93). 

One of the possible modifications oi M^^^^^^ is given by (4.81). 

Conjecture 3.18. All these estimates hold with 5 = 0. 

4 Superstrong magnetic field 
4.1 Preliminary analysis 

In this section we consider cases of very strong magnetic field^^^ 

(4.1) eo/?"^ < fi< Qh~^ 
and superstrong magnetic field 

(4.2) Co/)-' < /i 

In the latter case operator needs to be modified: we will replace V by 
V — fih^F where constant 3 will be specified later; so operator (2.1) is 
replaced by 

(4.3) ^= Yl Pjg^'^Pk-f^h^F+V, with Pj = hDj-fiVj 

l<J,k<2 

Then 



^^' Note that it is more narrow definition than (3.1) in the previous section. 
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(4.4) We need only to consider eigenvalues An(?7) of model operator L{rj) 
defined by (1.26) with n < N := A/(eo). 

Really, as n = 0, 1, 2, ... 

(4.5)d infAD,„(r7)>(2n + l) (n > 0) 

(4.5)^ infAN,n(77)>(2n-l), (n > 1). 

V 

due to proposition 6.1 and our operator is modelled by {^h)~^L{r]). 
However then 

(4.6) \X.,4r]) - X,,,{r])\ > e{N) Vm ^ n < A/ Vr/ 

with subscript * denoting either D or N and we can decompose into eigen- 
functions v^^nif^xi.^i) of the model operator. So basically we consider 
one-dimensional diagonal operator perturbed and we can diagonalize it to 
/i/?F(0, X2)An{x2, hD2) with 

(4.7) A(X2,/)D2) := (A,,,(nD2) -3 - (/i/))-^ 1/1/(0, X2)) + ... 
where in this section 

(4.8) h = ^i-h'2; 



preliminary we scaled xi 1— )■ h~^xi. Let us ignore perturbation in (4.7); later 
in subsection (4.5) we will show that it could be made very small, albeit 1/1/ 
will be modified. 

Remark 4-1- In (4.7) we replaced xi by 0. Actually the better approximation 
would come by replacing xi by fi^^hD2 but need in it will not come out 
instantly. 

Then we need to consider only scalar operators (4.7) and we are looking 
at fixed vicinity of some point X2. What we can say about .^2 apart of 
^2 > — Q (as X^n — > +00 as .^2 — ^ —00)? 

There are two possibilities: either we would be inside of the spectral gap 
if not a boundary, or not. 
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4.2 "Almost spectral gaps" 

Assume first that we would be inside of the spectral gap if not a boundary: 

(4.9) |(2m + l-3)/i/?F+ \/-r| >eo/i/? Vm = 0,1,2,... 

where |r| < C and usually r = 0. 

We know (see subsection 13.5.1 [Ivr2]) that inside of domain (4.9) means 
a spectral gap. However it is not the case near boundary: looking at behavior 
of \*,n{^2)fJ'h interpreted as Landau levels we can conjecture that only lower 
spectral gap (— oo, .) survives and even it shrinks for Neumann boundary 
condition. 

As X*_n{(,2) —^ (2n + l)fih as ^2 -^ +00 due to proposition 6.1 again and 
we conclude that 

(4.10) Under condition (4.9) operators in question are elliptic as ^2 > Q 
and we need to consider only a compact interval |^2| ^ Q- 

Note that as 

(4.11)„ Kn{^2)-d-{l^h)-'W>efih V6 

then operator An is elliptic for all .^2 and the contribution of it and interval 
(— cxD, r) to the asymptotics is negligible. In frames of (4.9) one can rewrite 
above condition as 

(4.12)d,„ {2n + l-^-e)i2hF+V-T>0, 

(4.12)^ „ {X*^^-^-e)fihF+V-T>0, 

for Dirichlet and Neumann boundary conditions respectively where 

(4.13) A*N „ = min An n{v) e (2m - 1, 2m + 1). 

One can then concludes easily that 

Theorem 4.2. Let one of assumptions (4.1), (4.2) be fulfilled. Further, 
let (4.12)pQ or (4.12)|\jq be fulfilled with r = (matching to the boundary 
condition). 
Then 

(4.14) |e(x,x,0)| < Cs/i"' 
with arbitrarily large exponent s. 
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Assume now that condition (4.9) is fulfilled for all m but (4.12)^ q fails. 
Then there is a large difference between Dirichlet and Neumann cases because 
only in the former case in virtue of proposition 6.2 

(4.15) a; „ < -e(A/, Co) Vn < A/ V^ : |6l < Q 

where here and below A'^ „ := 9^2^*," ^^'^- -^*^^ Neumann boundary condition 
we assume that 

(4.16) |A'n,,(6)I < eo ^ \{Xu,,{^2)-^)f^hF-W\ > eo^ih. 

So, we arrive to the remainder estimate 0(1); transition from Tauberian 
expression to magnetic Weyl one is trivial: 

Theorem 4.3. Let one of assumptions (4.1), (4.2) he fulfilled. Further, let 
(4.9) and in case of Neumann boundary (4.16) he fulfilled. Then R^^ = 0(1). 

However there are other mechanisms to break degeneracy for Neumann 
problem: first of all we know (see proposition 6.3) that 

(4.17) |A'n,„(6)| < eo ^ A';„(6) > eo Vn = 0, 1, 2, ... 

The first non-degeneracy assumption linked to derivatives with respect to 
X2 is 

(4.18) \{XN.n{^2)-i)l^h-W\ < eoi^h, |A'n,„(6)| < ^o ^ 

\d^,W\>eo Vm = 0,l,2, ..., 

where as usual 1/1/ = —VF~^; the next one is 

(4.19)^ \{\N.n{^2)-})f^h-W\<eofih, 

|An,„(6)I + I<9x, 1/1/1 <eo =^ ±dlW>eo Vn = 0,l,2,... 

Theorem 4.4. Let one of condition (4.1), (4.2) be fulfilled. Further, let 
conditions (2.60) and (4.9) be fulfilled in B{x, 1) and letip e ^^ he supported 
in e(x, |)n{x :xi < e'}. 

Then Rjj^* = 0(1) under one of the assumptions (4.18), (4.19) and 
pMW _ Q^i |og/j|)^ under assumption (4.19)~. 
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Proof. As /i/) X 1 we are in frames of subsection 5.1.3 [Ivr2] and arrive 
immediately to conclusion of theorem. 

Consider case (4.2) of superstrong magnetic field and adapt arguments 
of subsection 5.1.3 [Ivr2]. Then first we consider e(l, £)-admissible (with 
respect to (x2,.^2)) with p = edA'^^l partition of zone {|A'f^^| > C(/i/))^2|. 
Then for operator (4.7) 

(4.20) 7, X hp-\ T* X p-\ 

as propagation velocities with respect to X2,^2 are x p and 0((/i/?)~-^) 
respectively. Then contribution of an element to the Tauberian remainder 
does not exceed Chr^p x T^: x T* ^-^ ':< C and the total contribution of zone 
in question does not exceed C j p^^ dp^ C as operator is elliptic unless p 
has a specific magnitude. 

In the remaining zone {|A'f,j „| < C{fih)^2 let us introduce e(p, £)-admissible 
partition with 



1 



(4.21) p = {\X'^J' + {phr'\d.,Wfy + P. i={ph)-^P. P = /i-^. 
Then for operator (4.7) divided by ph 

(4.22) T,^hp-\ T*^{ph)K 

as propagation velocity with respect to {x2,{ph)2^2) is x p. Therefore 
contribution of (p, i) element to the Tauberian remainder does not exceed 
Ch~^p£ X T* X 7"*^-*^ X C and the total contribution does not exceed 
C J p~^ dp which is 0(1) under one of the assumption (4.18), (4.19) as in 
the former case p x pi = (p/))~2 and in the latter case operator is elliptic 
unless p has a fixed magnitude. Further, we get 0{\ log h\) under assumption 
(4.19)~ as integral is taken from p to pi. 

Transition from Tauberian expression to magnetic Weyl expression is 
trivial. D 

Remark 4-5. Following section 5.1.3 [Ivr2] we would be able to establish 
remainder estimate 0(1) even under assumption ((4.19) . 

Remark 4-6. (i) The principal part of asymptotics in question is 0{ph^^) 
and as 

(4.23) {! + }) + F-\V-T)<-e 
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it is X fih^^. 

(ii) As Neumann boundary condition is considered and 

(4.24) fih{l + 2,)F+V -T>efxh 

the principal part of asymptotics in question is 0{h~^) and as 

(4.25) l^hF{X*^o + d)+V-T< -ei^h 
, 11 

it IS X /l ■" = /i2/7 2. 

(iii) As 

(4.26) fih{2m-l + k)F+V -Ti<-efih, 

Hh{2m + 1 + k)F +V -T2>-enh, n < T2, 

then the increment of principal part of asymptotics calculated as r increases 
from Ti to T2 is 0{h~^) and is x h~-^ as \t2 — Ti| x fih. Recall that without 
boundary it would be 0{h^). 

Problem 4-7. Investigate under weaker non- degeneracy conditions than 

(4.19)^. 

4.3 Landau level 

4.3.1 Preliminary analysis 

Now consider case when condition (4.9) is broken. Without any loss of the 
generality we can assume that 

(4.27) Condition (4.9) breaks only for one value oi m = m. 

Really, if fih > C this assumption would be fulfilled automatically; otherwise 
we achieve it by considering e-admissible partition of unity. 

Further, without any loss of the generality we can assume that 

(4.28) 3 = 2m + l. 

Really, let i^ = \} — (2m + 1)| < 1. If ufih < C then term ±unhF could 
be absorbed into V; if ufih > C then dividing operator by ufih we find 
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ourselves in the frames of the previous subsection with /inew = ji^h^^u^^ 
and /inew = fi'^hiu'i < 1. 

Then for all n ^ m the previous arguments work and we need to consider 
n = m only. Further, as n = m these arguments work in zone {^2 < Q} and 
we need to consider only zone {Q < .^2 < eo/u}- 

As symbol operator fihAn with An defined by (4.7) has a either first or 
second derivative with respect to X2 disjoint from 0, but derivatives with 
respect to ^2 could be pretty small, we will see that we cannot skip without 
reservation 0{h) terms like those with an extra factor xi even as we derive 
Tauberian remainder estimate but we in our assumptions will be able in 
Tauberian arguments only to skip 0{h^) terms like those with an extra 
factor xl so we need to modify (4.7) a bit^^-*. To do this we first drag out 
F(x) out of [A — r) leaving inside 

(4.29) (/)Di)2 + {hD2 - fix.f - l/l/(xi, X2) + ■ ■ ■ = 

(/jDi)2 + {hD2 - fiXif - W{i2-^hD2,X2)- 

{d,, W){fi-HD2, X2)(xi - fi-HD2) + ■■■ = 

[hD^f + {hD2 - /iXi + ^/X-^(9,,l/l/)(/i-^/7D2,X2))' - W[fl-^hD2,X2) + ... 

with W = —{y — t)F~^ resulting after shift 

(4.30) hD2 ^ hD2 - fixi - -fi'\d^, W){fx-^hD2, X2) 
in 

(4.31) /i/?A(x2, hD2) = i2hX,,n{hD2) - W{h^D2,X2) + ... . 
There could also be 0(xi) terms^^^ 

(4.32) a(x2)xi (^h^Dl - (/7D2 - fix^f^ ; 

then replacing xi 1— )■ xi + /3(x2)x^ we also need to replace 

/?Di ^ (1 - 2/3x1 + 0{xf)) hDi and hD2 ^ {hD2 + 0{xf)hDi) 

^^) Actually effectively xi could be marginally more than h. 

^^^ After division by F intensity of magnetic field is 1 so there could be no terms (4.32) 
with "— " replaced by "+" . 
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so modulo O(x^) 
(/)Di)2 + {hD2 - fixif ^ (1 - 4/3xi)(/)Di)2 + (1 + 4/3xi)(/)D2 - ^xi)2+ 



4 

2 I /Ln I ^ a . .~1 U.2 ri2\2 



and choosing /3 = — |q; we arrive to 

(/)Di)2 + {hD2 - fixi + iPii'^h'Dl) 
and thus we need to replace (4.31) by 

(4.33) iih{An{x2, hD2) - 2n - l) = 

l2h{X,,,{hD2 + 2/3^^01) - 2n - 1) - 1/1/(^2 02, X2(l+) + 



Remark 4-8. In the case of superstrong magnetic field we need to consider 
only zone 

(4.34) {6, |A,,„-2n-l|< Qifih)-'} 
otherwise operator is elhptic. 

4.3.2 Non-degenerate case 

Consider first case 

(4.35) |(2n + 1 - 3)^/)F + \/| < eo ^ iVax^/F"^! > eo 

Vn = 0,l,2,... 

Then obviously one can take T* x h. Really, speed with respect to ^2 is 
x|9,,l/l/|xl. 

Meanwhile consider time direction in which ,^2 increases; we can select it 
due to assumption (4.35). To reach ^2 = e^~"^ we need time x h~^iih as for 
^2^ h~^ speed with respect to ^2 is {fih)~^. Consider evolution with respect 
to X2; speed with respect to X2 does not exceed C/i/7|A'^ „| + Ch. Then we 
can take for sure take 

(4.36) T* = r(6) X m\n{{f,h\Kj)-\h-') 
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Then contribution of B[x, 1) fl {x : xi < e'} to the remainder does not 
exceed 

(4.37) cy"r^^||jc/6 < J{h+Cf,h\Kj) c/6 

as integrand in the left-hand expression does not exceed Ch as |A'^„| < 
h{fj,h)~^ and it does not exceed C/i/) |A'^„| otherwise. As A*,n(.^2) is monotone 
(at least for ,^2 > C) and integral is taken over ,^2 < e/i^-'- also satisfying 
(4.34) the right-hand expression does not exceed C. 
Thus 

(4.38) Under conditions (4.1) or (4.2), (4.28) and (4.35) fulfilled in B{x, 1) 
contribution of 6(x, |)n{x : xi < e'}n{.^2 > C} to the Tauberian remainder 
does not exceed C. 



4.3.3 Generic case 

Now look what happens as (4.35) is also broken and replaced by conditions 

(4.39)^ |(2n + l-3)/i/jF+\/| + |VaxVF-^| <eo =^ 
and 



T'^lxVF-^>eo Vn = 0,1,2, 



(4.40)^ \{2n + l-2,)fihF+V\ + \VdxVF-^\<eo ^ 

^VVF~^>eo Vn = 0,1,2,... 

Inner zone. We define this zone preliminary by 

'Vinn = {6 > Co(log;U)2} 

which means not only that /i/)(A* „ — 2n — 1) and fihX'^ „ are negligible but that 
they remain so even if we replace .^2 by (1 — e).^2- Then under assumption 

(4.40) we can take 

(4.41) 7, xmin(e2"\nr2), 

(4.42) £x 11/14,1 
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Really, the propagation speed with respect to X2 is x /i but as scale with 
respect to ^2 is x ^2 there, uncertainty principle means that hT • ^2 > ^■ 
Similarly, propagation speed with respect to ^2 is x £, scale with respect to 
X2 is X £ and uncertainty principle means that iT ■ i > h. 

Consider propagation in the time direction sign(l/l/x2)t > in which ,^2 
increases. Then 



sign(l/14J^log|l/14J = -ni/14,sign |l/14J-^l/14, 



so log I l/l/xj also increases provided Wx^ and lVx2X2 have opposite signs. Then 
we can take 7"* = h~^. This would lead to contribution of X\„„ fl {,^2^ ^ ^} 
to the Tauberian remainder not exceeding 

c jj ^~'y^) "'^'"'''' = ^ // ""'"^^2^'' ^^"'^ "'^"'^^ - ^- 

On the other hand, in X\nn H {^2^ ^ ^} we can take 7"^, x ,^2^^ and 7"* x ^2 
and contribution of this subzone to the Tauberian remainder not exceeding 

C jj ^"'^lly di2dx2 = CJ i^^ d^2 < C. 

On the other hand, as (4.39) , (4.40) with the same signs are fulfilled 
we should chose between T* x ,^2 and T* x h^-^i. So, as £ > {hC,2)^ we 
can take T^ = M~^, T* = M~^ and contribution of this subzone to the 
Tauberian remainder does not exceed 

c jj h-'^^^d^2dx2 = c j r' di < c 

as in the case of the same signs £^ + ^2^ must have the same magnitude or 
operator is elliptic; so in this subzone i must have the same magnitude. 
As £ < (/i^2)2 we can take T* = ^^^, but the same ellipticity argument 

means that actually we can upgrade 7"* x /i to 7"* x /i 2^| contribution of 
this subzone to the Tauberian remainder does not exceed 

1 "^*(6) ,. , _ r f c-i 



C jj ^ y^ ^6c/x2 = C j ^,' d^2 < C 

as in this subzone i must have the same magnitude. 
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We need to extend X\nr, up to 

(4.43) A'inn := {6 > e7} 
with 

(4.44) ^f = (\og{fih'2)-{2n + 2)\og\og{fih'2)±C 

describing a zone where fJ'h\\'^,^{C,2)\ < e^ is dominated by \d^2^{^ih,X2)\. 
Then the same arguments work there albeit scale with respect to ^2 is now 
^2 — ^2"' ^2"' ^^^ defined by the same formula with C replaced by C/2. 

Further, these arguments remain valid for ,^2 ^ ■^2^ provided A'^ „ and 
— 1/1/xi have the same sign. 

Thus we arrive to 

Proposition 4.9. (i) Under conditions (4.1) or (4.2), (4.28) and (4.39)"^, 
(4.40) (with the same or opposite signs) fulfilled in B[x, 1) contribution of 
B{x, |) n A'inn n {.^2 ^ C"} to thc Taubcrian remainder does not exceed C . 

(a) Further as A'^ „ and — Wx^ have the same signs^^^ (i) remains true for 

'^inn = {^2 > 'C2 } ■ 

Boundary zone. Consider now boundary zone temporarily introduced as 

'Abound '■= {^2 < ^'^2} 

which actually should be intersected with zone (4.34). 
Then we can take 

(4.45) T,^m\uh{i2[tih\Kjy\r^)- 

Really, in this zone propagation speed with respect to X2 is x /i/?|A'^ „| and 
scale with respect to ^2 is x ^2^ to keep A'^ „ of the fixed magnitude. 

As £ > (^2^"^Ai/'|A'^ „|) ^ consider propagation in the time direction of ^2 
increasing. Then 

(4.46) 7*xmin((^/)|A:,„|)-\ri£) 

^^ Which means that under Dirichlet or Neumann boundary condition (4.40) and 
(4.40)^ respectively must be fulfiUcd. 
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provided | l/l/x2 1 also increases in the same direction because speed with 
respect to X2 is x (/i/j|A'^ „|) as long as we remain in A'bound but then drops 
to 0{h) outside of A'bound but X2 can go in the opposite direction there. Our 
extra assumption means exactly that 

(4.47) Wx2X2 ^^(^ K,n have the same signs^^^. 

Then contribution of the subzone in question to the Tauberian remainder 
does not exceed 

X C J {{f,h\Kj^2Y^ + f^2{f^h\Kjr') di2 X 

c{i^h\KM2)k2\^^, + ch[^^h\K^^)-\^-^^ X c 

where ^2 is defined from (4.34) as solution to yu/?|A*,n — (2n + 1)| = C. 
Consider now subzone 

where the second condition is due to uncertainty principle and scaling with 
respect to .^2- 

Then 7"* is defined by (4.45) with (i replaced by [^2^iih\\'^ „|) ^ as we can 
go into direction of increasing ^2- Then contribution of this subzone to the 
Tauberian remainder does not exceed 

^//^ '^H) '^^'"'''' ^ Chjdi2 « c 

again. Further, in subzone {£ < ^1.^2} one can take T* '^ ^2^ and its 
contribution is ^ 1 as well. 
Consider now case when 



(4.48) H/x2X2 and A'^ „ have the opposite signs 



16) 



15) Qj. equivalently that under Dirichlet or Neumann boundary condition (2.39)" and 
(2.39) respectively must be fulfilled. 

16) Qj. equivalently that under Dirichlet or Neumann boundary condition (2.39) and 
(2.39) respectively must be fulfilled. 
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Then as £ > £* we can do nothing better than 

(4.49) T* = i{f^h\KJ)-\ 

Then contribution of the subzone in question to the Tauberian remainder 
does not exceed 

X C / c/^ X Clog/i. 

On the other hand, for i < i* we can improve (4.49) to 

(4.50) T*=e~'{f^h\KJ)-'-'^'. 

Then contribution of the subzone in question to the Tauberian remainder 
does not exceed 

X C / c/5 X Clog/i. 
FinaUy, extending A'bound to 

(4.51) A'bound = {6<e2"} 

comes with no cost at all as scale with respect to .^2 is ^2^^ x (,^2^' — ,^2) 

anyway. 

So we arrive to 

Proposition 4.10. (i) Under conditions (4.1) or (4.2), (4.28) and (4.39)"^, 
(4.40) (with the same or opposite signs) fulfilled in B{x, 1) contribution 
of B{x, |) n A:'bound n {^2 > C} to the Tauberian remainder does not exceed 
Clog/i. 

(ii) Further as \'^ ^ and l/l/x2X2 have the same signs^^^ this contribution does 
not exceed C . 

As A'^ „ and — W^^ have the same signs^^-* we are done. Otherwise we 
need to consider 
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Transitional zone. A'trans is defined by 

(4.52) e^ < 6 < e7 

with ^2 defined by (4.44) and in this zone /i/jA'^ „ has magnitude h and 
///jA'^' „ has magnitude /i^2 ^ ^^ and we conclude that 

Proposition 4.11. Under conditions (4.1) or (4.2), (4.28) anc? (4.39) , 
(4.40) (with the same or opposite signs) fulfilled in B{x, 1) contribution 
of B(x, I) n A'trans H {^2 > C} to t/ic Tttubcrian remainder does not exceed 
Clog/i. 

(a) Further as \'l ^ and 1^x2x2 have the same signs^^^ this contribution does 
not exceed C. 

Synthesis So, we proved that the total Tauberian remainder does not 
exceed Clog/i. To upgrade it to C we must assume that l/l/x2X2; A'^ „ and 
— 1/1/xi have the same signs as A'^ „ and \'l „ have opposite signs. So we proved 



Proposition 4.12. Under conditions (4.1) or (4.2), (4.28) and (4.39) 



± 



Vlt 



(4.40) (with the same or opposite signs) fulfilled in B{x, 1) contribution of 
B(x, |) n {.^2 > C} to the Tauberian remainder does not exceed C log//. 

(ii) Further as l/l/x2X2; A'^ „ and —W^^ have the same signs^^^ this contribution 
does not exceed C . 



4.4 From Tauberian to magnetic Weyl 

In this subsection T = T^^jJ but we will calculate with (5 = because if we 
consider time interval [T, 2T] with 7" > T^, in the estimate of its contribution 
there will be an extra factor 7~^^ 7"^ and the summation with respect to T 
will result in the same answer albeit with T = T^. 

We apply the successive approximation method to operator (4.7) 

4.4.1 Inner zone 

Replacing X2 by 3/2 leads to a contribution of (i, p) element to an error does 
not exceeding 

(4.53) Ch-^pi X — ^ = Ch-^piTl 

It 
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as propagation speed with respect to X2 is 0{h) in the time scale used here 
and we use time scale compatible with the choice of T* in the previous 
subsection. 

Under non-degeneracy condition (2.8) we can take T^, = h resulting in 
the contribution of A'jnn equal to 0(1) as p x h^-^. 

In more general generic case we estimate a contribution of [i, p) element 
to an error by 

ft/ 72 

(4.54) Ch-^pi X — — = C^2f Tl X Ch-^i2f min (^2^^ h^r'^) 

It 

where we plugged T^, defined by (4.41); an extra factor d, appears as we 
replace X2 by 3/2 because 8x2 W = 0(£) and at this moment we take original 
"narrow" X\„n with scale /) x ^2 in ^2- 

Therefore the contribution of A'jnn to an error does not exceed 

(4.55) Ch-^ ffm\n{^f,h^r'^)id£d^2 x C [^2^ d^2 x Clog/i. 

In the remaining part of A'jnn defined by (4.44) we must replace ,^2 by 
p X ^2 ~ ^2 as the scaling is concerned and it redefines correspondingly 7"^, 
and instead of right-hand expression in (4.55) we get C J{^2 — ^2)"^ d^2 
which is also x C log p. 

4.4.2 Boundary zone 

Due to the standard arguments contribution of subzone {^2 < Q} to the 
error is either 0{\ogp) or 0(1) depending if there is a saddle point or no; 
one can see easily that this saddle points cannot come from AN,m but only 
from A|\i,n with n ^ m and it is possible only as Xh x 1. 

Consider subzone {^2 > Oq}. Here we need to replace (4.53) by 

(4.56) Ch-'pi X ^^M^^^Jll = Ch-'ph\KJpiT^ 

as propagation speed with respect to X2 is x /i/)|A'^^. Plugging under 
condition (2.8) i = 1, p = ^. 



-1 
2 



r, = nmin(6(/i/)|A:,„|)-\l) 

we arrive to 

Cph\Kjpmin{e2{ph\Kj)-\l) 
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and contribution of both subzones {.^2 > C, f^h\\'^ ^\ > ^2} and {Ch < 
l^h\KJ < 6^} are 0(1). 

In the same manner (4.54) is replaced by 

(4.57) Ch-'pf X ^l^^^yi! X Ch-'f,h\KJpfT! 

and plugging 

we arrive to 

C^h\KJpfvr^\n{e2iph\Kj)-\r') 

and summation with respect to i results in 

CphiKjpt'-^ = Cpi2 

as £* = ^2^ ph\\'^,j^\ and summation with respect to .^2 results in C^2^ ^ 
C\ogp. 

4.4.3 Transitional zone 

In this zone we can return to (4.53) and plug T^ = h under condition (2.8) 
resulting in 0(1). 

Further we can return to (4.54) and plug in the generic case either 
T^ = min(M~^,.^2^"'") or 7"^, = min(M~^, |^2~^2|~''") resulting after summation 
with respect to £ in C^2^p or C\^2 — ^2\^^P respectively. In the former case 
we instantly get 0(1) while in the latter after summation with respect to 
subzone {|'C2 — 'C2I > ^} we arrive to O(log;u) while contribution of subzone 
{IC2 — 'C2I < ^. ^ < ^} to asymptotics is 0(1). 

4.5 Justification: reduction to model 
operator 

4.5.1 Reduction: step 1 

The problem however is that the "kinetic" part of our operator is not exactly 
h^Dl + {hD2 — pxi)'^ but is different. We need to improve construction of 
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subsection 2.1). First we can assume that Vi = and that g-''^ = 5jk, V2 = 0, 
t)i V2 = 1 as xi = 0. Further, assume that 

(4.58) g'^ = 0{x^). V2 = x, + 0{x!^+') 

with k > 1. Then changing xi 1— )■ xi, X2 h- )■ X2 + a{x2)xi^^ we can achieve 
g^^ = 0{x-^^^), simultaneously preserving the second relation. Then re- 
defining Xi := V2(x) preserves the first relation (with k + 1 instead of k). 
Continuing we can make k as large as we wish and thus k = 00. 

One needs to remember that g^-^ and g^^ differ from 1 by 0(xi). However 
as we consider operator F~^A with the intensity 1 of magnetic field we 
conclude that g^^g^^ = 1 + 0{x^). 

Our goal is to get rid off this perturbation which we rewrite as 

(4.59) P™ := /i^ (n^DiaDi - a{nD2 - x{)a{nD2 - x{) 
Consider Poisson brackets 

(4.60) l{(^2^(^2-xi)2),a + /3ei} 
and note that for coefficients at ^1 equal i.e. 

(4.61) a,, - (xi - 6)/3x2 = 
it is equal to 

(4.62) I3^^^l_^i3 + a,;)i^y,^-i2)- 

To make it equal to (4.59) modulo ^l + (xi — ^2)^ one needs to satisfy 

(4.63) a,, = 2a(xi - 6) - P.A>^i - 6) - /3- 
Compatibility with (4.61) requires 

(4.64) {dl + 9^J ((xi - 6)/3) = 2d., ((xi - 6)^) . 

Consider instead it without d.^i with an extra condition "/3 = as xi = 0'' 
We can solve this one-dimensional equation 

(4.65) dl ((xi - 6)/3) = 2d., ((xi - 6)^) , 
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then 

(4.66) /3(xi, .) = 2(xi - 6)-' r W - 6)^(x(, .)dx[ + p 

with p which does not depend on Xi; extra condition "/3 = as xi = 0" 
means that 

(4.67) P = -^2' [ \x'i - ^2Hx[, .)dx[ 

Jo 

which gives us smooth function /3 and a. 

Remark 4-13. As ax26 = as xi = ,^2 = due to (4.63) and properties of 
a, /3 we can select a such that a^2 — there as well. This would enable 
further calculations. 

Then 

(4.68) ^-{{e, + {^2-x,Y).a + ^^,} = 

(/3x, - ct) {e, + (xi - 6)') + ^(^1 - {xi - 6)') . 

Let us pass to operators; to do this we consider corresponding /i-quantizations, 
divide by —iH and multiply by /i^; then 

(4.69) -^[A,C'"]=V'"-r^, 
with 

(4.70) A:=h^Dl + [fiD2-xif, 

(4.71) £:=i(« + /3^,), 

(4.72) Pi = (a - /3.,a) {il + (xi - 6)') 

where as we consider Weyl quantizations and symmetrized products after 
multiplication by /i^ the error is 0{^^1t?) = 0{h^). 

One can see easily that in our settings (T = /3 = 0asxi = and (3, a^^ 
has of the second order as xi = ,^1 = ,^2 = while a^^ has of the third 
order as xi = ^1 = ,^2 = 0. Also note that due to (4.64) /3xi — o" = as 
xi = ^2 so 
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(4.73) symbol (a — /^xj is divisible by (xi — ^2)- 
Consider transformation of A + V*' by exp(//i^-^£™): 

(4.74) exp(/r '£"') {A + P") exp(-/r Y™) = 

A + v" + ^[A + v'', r] + ^(^)'[[/^ + v", rir] + ... 

/n 2 /ft 

where as one can see easily we left out terms which are o(/)^) and due to 

(4.69) 

(4.75) =A + V'^ + \i^)[V'",r] + \{^)[V'^,C'\ 

2 in 2 ifi 

Consider first the symbol of the third term 

- ^-{a, a + /36}(ei - (xi - 6)') - \<y{ei - (xi - i2)\ a + /J^i}; 
as 

^{ei-(xi-6)',a+/36} = K+/3x,(xi-6))6+/3xiei+(xi-6)(«x,+/3) = 
2/3x,(xi - 6)^1 + ^(^1 + (xi - 6)') - (^ - /3xJ(ei - (xi - 6)') 
we conclude that 

-J{P,a + /3ei}=(-^{^,« + /36} + ^^(^-/3xj)(ei-(xi-6)')- 

a/3x,(xi - 6)6 - \^\ei + (xi - 6)'). 
Consider now the symbol of the fourth term in (4.75): 

-J{Pi,a + /36} = 

-^{(^-/3xJ,«+/36}(ei+(xi-6)')-J(^-/3xJ{ei+(xi-6)',«+/36}; 
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as we know the last Poisson bracket we conclude that 
-Jm,« + /36}=(-J{(^-/3xJ,« + /3ei} + ^(^-/3xJ')(ei + (xi-6f) 

Then the sum of symbols of the third and the fourth terms in (4.75) is equal 
to 

(4.76) - ^{a, a + /^CiK^i - (^i - 6)') - ^/3x,(xi - 6)6+ 

{-\{{^ - /3xJ, a + /36} + \{-^' + (^ - /3xJ')) (Ci + {xi - 6)'). 

Remark 4-^4- However correction (4.59) is not a correct one as we need 
^11^22 ^ -^ g^j^j ^ more precisely is Vi + \cP-A as g^^ = 1 + a + \cp- + ... , 

We can skip terms here which has of degree 5 at xi = 6 = 6 = 0? this 
takes care of the middle term. 

By the same reason we can skip in 

--{a, a + /36} = ^f^xjtt^ + 4^x2/366 + 4^xi/3 

middle term where we used that a = (j{x). 

Consider two remaining terms (which have of the second order at 
xi = 6 = and also equal as xi = ) 

(4.77) ai := -ay^^a^.^ + -a^,l3; 

then we can repeat the same procedure as before thus replacing the whole 
first term in (4.76) by 

(^i-/3ixJ(ei + (^i-6)') 

with Pi matching o"i. One can see easily that then neither Vi nor the 
last term in (4.76) would change modulo 0{h^) and therefore we arrive to 
operator A + P™ with 

(4.78) V2 = 

((a-/3.J + ai-/3i.,-^{(a-/3.J,a+/36} + ^(^-/3xJ')(ei+(^i-6)'). 
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In this factor the first term has of the first order, other terms have of 
the second order or higher and we can skip those which are higher. We can 
also skip terms divisible by (xi — ^2)^ as those terms are 0{H) on energy 
levels in question leading to 0{h^) error in the final answer. Due to (4.73) 
we arrive to 

{a - /3xJ + (^1 - /3ixJ - J{(^ - /3xJ, a + f3^i} 

where the last term is 

- J{(^ - /3xJ, a} -\{{a- /3xJ, ^i + (^ " Mx^P■ 

One can prove easily that due to (4.63), (4.73) (which is valid for (Xi, Pi as 
well) the remaining terms vanish as xi = ^2- 

Therefore as we skip terms containing factor (xi — .^2)^, what is left of 
perturbation is V^ with 

(4.79) V3 := p(x2, 6)(xi - 6) (^1 + (xi - 6)') , 

so, V^ is a symmetric product of p™, (xi — ^^^hD2) and A (as factor fi^ 
should be remembered). 

4.5.2 Reduction: step 2 

This perturbation is of the same magnitude as the original perturbation V^ 
(albeit h~^ in estimates does not pop-up) but it is much easier to handle. 
If we had no boundary^^^ we would use C = uj{x2, ^2)^1(^1 + {xi — ^lY) to 
eliminate it generating lower order terms but closer to factor ^1 forbids it 
unless accompanied by factor xi. Still replacing p(x2,^2) in (4.79) by 

p{x2, 6) = p{x2, 0) - p^^xi + p^2(xi - 6) + + ... 

we can eliminate the second term (producing exactly such component in V3) 
and skip the third and all terms as producing 0{h^) contribution. So 

Remark 4-15. Without any loss of the generality one can assume that p = 
p(x2) in (4.79). 



17) Or were on the distance at least h^ ^ from it. 
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Still it does not help us much and we cannot use obvious transformation 
by exp{iH^^C") with C = uj{x2)[^l + (xi — ^2)^) due to trouble in the next 
commutators. Instead we recall that basically our operator (4.70) multiplied 
by /i^ should be equal W + fih^ with constant 3. Actually one needs to insert 
there also —hDt but at energy level it is really small and we can actually 
by multiplication by (/ + (?) make our evolution equation linear with respect 
to hDf. We leave absolutely standard detailed arguments to the reader. 

So now perturbation becomes 

(4.80) r: = fihp{x2){xi - nD2)}> + p(x2)(xi - nD2)i 

and it is completely different game as we can get rid of the first term using 
transformation by exp{ih~^ C") with C = fj,hu{x2)2i as one can see easily 
that there will be no problems with the next commutators. Remaining 
perturbation is of the type (l/l/i(x, ^)(xi — ,^2)) and is of magnitude h. So 
we gained factor [nh)~^ which is important only in the case of superstrong 
magnetic field so we can skip this step for a very strong magnetic field. 

4.5.3 Reduction: step 3 

Note that the only part of new perturbation which is not necessarily 0{h^) 
comes from {Wi{x2,C,2){xi — ^2)) and we can use the same approach as 
before to eliminate it. 

Actually we can continue further. We can keep power of HDi below 
2 replacing fi^D^ by — (xi — fi^D2)^ + ^3 — fi^^W. We can eliminate HDi 
transforming by exp{iH^^C^) with C = p(xi — ^2)- Similarly, we can eliminate 
x^(xi — .^2)'' with p > q transforming by similar operator with each ^1 having 
cofactor Xi. Therefore having x^ with n > 2 we can decompose it into 
sum of terms of the type x^{xi — tiD2)'^{tlD2)"^'^^^ where either p > q or 
p = 0,q = lorp = 0,q = and eliminate the first type by the method 
we just discussed and the second one transforming by exp{iH~^C^) with 
C = C{x2, ^2)- It eliminates all powers - including those we skipped before 
as lesser than 0{h^). 

Remark 4-16. (i) So in the end we arrive to the model operator. Sure 1/1/ 
is going to be perturbed by 0{h^) albeit it may be because of smallness of 
Ci] so actually we get in the end l/l/eff(6,^2) = -( V/F)(^2,X2) + 0(^| + h^) 
albeit this does not change it basic properties. 
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(ii) Recall that in this subsection ^2 in contrast to the rest of section 
corresponds to fiD2 rather than hD2 and h = f}^. 

Conjecture 4.17. l/l/eff(6.^2) = -(^//^)(6.^2) + 0{i^ + H^) as ^ < h^ . 

4.6 Final results 

Thus we arrive to the final contribution of zone {,^2 < h~^} (or equivalently 
{xi < h'i^^}) to the asymptotics: 

5^(27rn)-V //e(/i/)(A,,„(6)-3)-M/eff(^6,X2))^(6r\x2)C(6^') a'X2C/6 

where transformation of ip is obvious. After change of variables xi = h^2 

^(27r/))-l /"/"e(/i/)(X,„(xir^) -3) - l/l/eff(xi,X2))V^(x)C(xi5-l) C/X1C/X2 

„\n 'J ■J 



n>0 



with e = n}-^ and h = {i^-^h)^ Here C e "^^^{[-1, 1]) equals 1 on [-|, §]. 
On the other hand, the final contribution of zone {^2 > h~^} (or equiva- 
lently {xi > h2~^}) to the asymptotics is 

5^(27r/7)-V II Q{lih{2n+l)-i) - l/l/eff(xi,X2))V'(x)(l - C(xi£-')) c/xic/x2 

n>0 "^"^ 

and the final answer is 

J](27r/7)-v/"/"e(/i/7(2n + l)-3) - l/l/eff(xi,X2))VWcyxicyx2+ 

n>0 "^"^ 

J](27r/7)-V //(e(^/)(A,,„(Xirl) -3) - l/l/eff(Xi,X2))- 

e(/i/7(2n + l-3) - l/l/eff(Xi,X2)))^(x)C(xi£-l)cyXiCyX2. 



n>0 



where now obviously we can take e = Ch\ log h\ and even e = Ch as we are 
in the spectral gap situation. 
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This matches (3.17) with 



(4.81) AT^^ 



bound 



5];(27r)-V [[{e{fih{K,n{x,h-')-^) - l/l/eff(xi,X2))- 
n -^ -^ 

e(/i/)(2n + l-3) - l/l/eff(xi,X2)))^(x)C(xi£-^)c/XiC/X2. 



Therefore we arrive to our final theorem 

Theorem 4.18. Letip G ^°^(X) he a fixed function with a compact support 
contained in the small vicinity of dX and let conditions (2.60), (2.7) and 
(2.41) he fulfilled there. Then 

(i) R^^ withU^^^^ defined hy (4.81) does not exceed Clog^; 

(a) This estimate could he improved to 0(1) in each of the following cases: 

(a) Assumption (2.8) is fulfilled; 

(h) Assumption (4.9) and in case of Neumann boundary problem one of 
conditions (4.16), (4.19) is fulfilled; 

(c) Assumption (4.9) fails hut there is Dirichlet boundary problem and both 
conditions (2.41)~ and (2.7)^ are fulfilled; 

(d) Assumption (4.9) fails but there is Neumann boundary problem and 
both conditions (2.41) and (2.7) are fulfilled. 

Remark 4-19. Here we redefined J^^^^^^^ including in it in contrast to what 
we used before (see f.e. theorem 2.6) W[xi,X2) rather than 1/1/(0, X2). 

We leave to the reader to prove that replacing H/(xi, X2) by 1/1/(0, X2) leads 
to an error 0(1) under condition (4.8) and 0{h^2\ \ogh\'-) under condition 
(2.41). 
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5 Generalizations 

5.1 Robin boundary value problem 

Consider boundary condition 

(5.1) {iJ2^Jg''Pk-cy)uU^ = 

J.k 

with real- valued a = a{x) > where z/ = (z^i, ... , z/^) is an inner normal. 

As a = we get Neumann boundary condition and in a formal limit 
a — )■ +00 we get Dirichlet boundary condition. 

Weak magnetic field Obviously arguments of section 2 remain valid. 
Recall that in that section we have not distinguished between Dirichlet and 
Neumann boundary conditions. 

Strong magnetic field Obviously arguments of section 3 remain valid. 
Further, under assumption a > Cq propagation theorem 3.15 remains valid; 
really then boundary problem satisfies Lopatinski condition. Therefore 
under this assumption one can derive the same remainder estimates as under 
Dirichlet boundary condition. 

Very strong and superstrong magnetic field On the contrary, in this 
case we are practically in frames of Neumann boundary condition: see 
remark 6.8. 

5.2 Boundary meets degeneration of V 

The following problem seems to be rather easy and straightforward: 

Problem 5.1. Get rid off condition (2.12) V < — e as ///) < 1 (we have not 
imposed this condition as fih > 1). 

One should apply rescaling method and follow arguments of subsec- 
tion 13.7.1 [Ivr2]. 
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5.3 Domains with corners 

The following problem is far more interesting and challenging and I strongly 
believe that the detailed analysis merits a publication: 

Problem 5.2. Consider domains with corners 7^ 0, vr, 27r assuming that coeffi- 
cients belong to "^"^(X), all other assumptions are fulfilled and in the corners 
non- degeneracy condition |Vaxl^| ^ 1 holds for each of two directions along 
dX. 

Without any loss of the generality one can assume that locally either 
X = {xi > 0, X2 > 0} or X = M^ \ {xi < 0, X2 < 0}. Then non- degeneracy 
condition means that 

(5.2) |l/l/xjxl, |l/l/xjxl. 
We also assume that (2.f2) is fulfilled. 

5.3.1 Weak magnetic field case 

One can use rescaling method with I = e|x| as long as /leff = hi~-^ < 1, 
A*efF = ni > I. In the inner zone X\„„ (defined exactly as before) magnetic 
drift is the only way to break periodicity and one can see easily that (3.21) 
should be replaced by 

(5.3) T, := h^-^r^ 

Then we are in frames of the weak magnetic field approach as long as 

(5.4) i>I:=max{Ci2-\i2h^'^). 

Obviously we can take T* = efi'j with 7 := min(|xi|, |x2|). To increase it 
consider dynamics on Figure 9. 

Selecting a proper direction (away from the nearest boundary) we can 
take there 7"* = efii and the contribution of A'lnn intersected with ^-element to 
the remainder does not exceed Ch^^l^ x {fii)~'^ x Cfi~^h~-^i and summation 
over all inner zone returns 0{fi^^h^^). 

Meanwhile contribution of ^-element intersected with the boundary zone 
'Abound and transitional zone A'trans does not exceed Cfi~-^h~^ixi~-^ as 7"* = e£ 
here; then summation with respect to these zones returns 0{fi~^h~^\ \ogh\). 
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(a) W^, > 0, 1A4, > (b) 144^ > 0, W^, < (c) W^, < 0, W>,, > (d) W^, < 0, 144, < 




(e) W,, < 0, W,, < (f) I44i > 0, 144, > (g) 144^ < 0, 144, > (h) 144, > 0, 144, < 

Figure 9: Dynamics in the corner with angle < vr (a)-(d) and > vr (e)-(h). 
Bold lines show hop movement and thin lines show drift movement which is 
along level lines of 1/1/. Consistent with Figure 3. 



One can get rid off logarithm noting that one can upgrade T* to T* = 
rr\\n{v~^i, 1) as we are trying to move away from the corner; recall that in 
the transitional zone A'trans speed in the opposite direction may be at most 
V := C{iih)ih~^ + C^~^ due to analysis of section 5.3[Ivr2] and /x£ > 1. 
Further, in the transitional zone we have an extra factor C{nh)^h~^ due to 
analysis of section 5.2 [Ivr2]. 

We are left with the corner zone Afcomer ■={><■ ^(x) < £} and in virtue 
of simple rescaling x i— )■ iix, /x i— ;■ 1 and /? i-> /i/i < 1 we can estimate its 
contribution to the remainder by 



(5.5) 



c^h-^F = Cfi-H-^ (1 + /Z?^"'') ; 



So, the total remainder estimate is also given by (5.5) which is 0{fi ^h ^) 
for 



(5.6) 



/i< /? 
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5.3.2 Strong magnetic field case 

To improve estimate (5.5) (our ultimate goal is 0{fi^^h~^) as usual) as (5.6) 
breaks one needs to invoke arguments of sections 2 and 3 and take into 
account dynamics near the corner (see Figure 9) in the same way we took 
into account dynamics near critical point of 1/1/ 1 ax (see Figure 4) before. 

Tauberian remainder estimate. As we mentioned in X\„n one can take 
T* X fii and 7"* = h^~^i~^. Therefore contribution of ^-element intersected 
with A'inn to the Tauberian remainder does not exceed 

(5.7) Ch-^f{l + fiT,)T*-^^ Cn-^h-\i + nh) 

where due to arguments of section 2 here and below in the estimates we can 
take (5 = in the definition of 7"*. Then summation over X\„„ fl {x, £(x) > €\ 
results in 0(^fi~^h~^ + \ \og h\y Moreover, in cases (b),(c), (e)-(h) one can 
take T* ^ fi and the result wiU be 0{fi^^h^^). 

The weak magnetic field arguments remain valid in A'bound and its contri- 
bution to the Tauberian remainder is 0{fi~^h^^). 

Meanwhile exactly as in section 2 contribution of ^-element intersected 
with A'trans to the Tauberian remainder does not exceed 

as T* X m\n[l, v~^i) where v is the propagation speed in the direction 
opposite to the hops and [jih)ih~^ is the width of Atrans; then the total 
contribution of A'trans H {x, £(x) > £} does not exceed 

Cfi-\^ih)h-^-^ {l + fih+v + fihvr^) ; 

meanwhile the contribution of A'trans H {x, Ch < i < i} does not exceed 
C{jih)^h~^~^(i as we take T* x ji"^ here where here and below 

(5.8) n:= max(/i-\/?^-'^). 

Let us find d, from T* = T^ but keep it greater than h i.e. 

(5.9) ^= Cmax(n,(i//j)5); 
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so the total contribution of A'trans H {x, £(x) > h} does not exceed 

Plugging V = [fih)3h~^ we arrive to 

(5.10) Cfi'H-^ + Ch-^ + Ch-'2-\fih)l 

where the last term is not greater than the first two as /i < /?~it. 

On the other hand, under Dirichlet boundary condition^*^^ we can take 
V = fi~-^ but then fii > 1 and we can take T* = 1 and we arrive to 



Remark 5.3. As it takes time (jj,h)3h^ to punch through A'trans and shift 
with respect to X2 will be v[jih)ih~^ we can improve T* as v[fj,h)3h~^ < i 
which is not needed in the case of the Dirichlet^^^ problem and does not 
help in the case of the Neumann problem. 

We are left with .^corner but its contribution does not exceed CiJih~^(? = 
Ch-^. 

Therefore remainder estimate is given by (5.10); if on the both sides of 
the corner Dirichlet^®-* ^^^ is given remainder estimate is 0{(i^^h^^ + h^^). 

Magnetic Weyl remainder estimate. Transition from Tauberian to 

magnetic Weyl formula is obvious in A'inn- To run successive approximations 
one needs to assume that 

(5.11) {fi~^ + vT,)T,<h'+^. 

This is definitely the case in A'bound and we need to consider A'trans only. 
Plugging 7"^, = hi~^ we arrive to i > ii with ii defined by 

(5.12) ii = C max{i2-^h-\ (i//?)^) 

which brings no extra terms in comparison with the Tauberian remainder. 
So in contrast to section 3 transition from Tauberian to magnetic Weyl 
formula does not increase an error. 



^^' And Robin boundary condition with a > e. 
19) 'We can obviously mix conditions. 
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5.3.3 Very strong magnetic field case 

We need also to consider case of very strong magnetic field case h^~^ < 
fj, < eh^^ under Dirichlet and improve remainder estimate 0{h^^) we got 
here. In this case analysis follows basically from the same arguments 
as in subsection 3.1. Then in formula (5.7) one should replace (/i/?)?/?"^ 
by 1 resulting in the remainder estimate 0{fi~^h~^ + \ hgfil'^). In this 
approach actually corner zone becomes A'comer = {><• ^{x) < ^} with i = 
Cfi'^ + Ch\ \ogh\'^/^) or even 1= C/j,'^. 

Problem 5.4- (i) Investigate very strong magnetic field case under Dirich- 
let^^-* and recover remainder estimate 0{fi^^h^^ + | \ogfi\'^); 

(ii) Improve it under Dirichlet boundary condition to 0{fi~^ h^^); is it 
possible? 

(iii) Could be it done under Robin boundary condition? 

5.3.4 Superstrong magnetic field case 

The similar arguments should work in the superstrong magnetic field case 
fi > eh~^ albeit we need to take i = fx~2+^h2, i = n~2h2\ \ogfi\'^^^ and 
i = C/i~2/)2 to recover remainder estimates 0{fi^), 0{\ log/i|'^) and 0(1) 
respectively and we arrive to 

Problem 5.5. (i) Investigate very strong magnetic field case under Dirich- 
let^^-* and recover remainder estimate 0(/i^); 

(ii) Improve it under Dirichlet boundary condition to 0{fi~^h~^ + | log/i|'^) 
and even to 0{fi^^h^^); is it possible? 

(iii) Could be it done under Robin boundary condition? 

5.4 Boundary meets field degeneration of F 

5.4.1 Set-up 

Even more challenging (and in my opinion worth both publication and 
degree) seems to be a problem 

Problem 5.6. Get rid off condition \F\ > e. 
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Our goal (not necessary achievable in all cases) is to derive asymptotics as 
sharp as if there was no boundary. We will assume that F x 'y(xy~^ where 
7 e ^°°, 7W = ^ |V7(x)| X 1, X G 9X,7(x) = ^ |Vex7(x)| x 1, 
V <0. 

Then we will have Xreg = {x : |7(x)| > 7 := Cfi^^''^} in the case of 
H < eh~'^ and degeneration zone Xjeg = {|7(x)| < 7} and then we will have 
also four subzones: 

(5.13) X,eg,inn = {x i 7(x) > 1, ^(x) != dist(x, SX) > C fl-'^ixy-"} , 

(5.14) X,eg,bound = {x : 7(x) > 7, ^(x) < Cfi-'^ixy-'-}, 

(5.15) Xdeg.inn = {x : 7(x) < 7, ^(x) > Cfx-'^'"" = C7}, 

(5.16) Xcorner = {x : 7(x) < 7, ^(x) < Q} 

as on Figure 10(a) and each zone should be analyzed separately. 

Remark 5.7. (i) Actually we need slightly increase 7 as /i is close to h~'^ 
and we need increase "slightly" the threshold for £(x) as 7(x) is close to 
(/i/?)-!/^'^-!) for fi > h^-\ 

(ii) When we consider modified Schrodinger operators with fih > 1 so 
that domain {x, 7(x) > C{fj.h)^^'^'^~^'} is no more forbidden, according 
to section 4 we need to change threshold for i{x) to "slightly more" than 

^-l/2;,l/2^(^)(l-.)/2_ 

(iii) When we consider modified Schrodinger operators with ^h'^ > 1 accord- 
ing to section 14.8 [Ivr2] we need to change threshold for 7 to "slightly more" 
than /x-i/2!./,i/2_ 

Without any loss of the generality one can assume that X = {xi > 0} 
and 51 = {x2 = 0} and F12 < as X2 > 0. Then hops here are to the left (x2 
decays) and drift down (xi decays) because dx2{—V)/Fi2 > 0. 

Consider X2 < 0; note that F12O changes sign and c?x2(~^)//"i2 does not if 
1/ is even and then hops are to the right and drift is down there (Figure 10(b)). 
On the other hand, F12O does not change sign and dx2{—V)/Fi2 changes if z/ 
is odd and then hops are to the left and drift is up there (Figure 10(c)). 
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(a) 



(b) V is even (c) i' is odd 



Figure 10: Zones (boundary zone is gray, degeneration zone is dotted) 
and dynamics near degeneration line of F (dashed). As U < we consider 



7' 



u-l 



with smooth 7 such that IV7I x 1. Consistent with Figure 3. 



5.4.2 Weak magnetic field case 

Using methods of the previous chapter and considering time direction in 
which drift up ne can prove easily that contribution of Xregjnn to the remainder 
is Oin'^/^h-^) at least as /^ < h^'". 

Furthermore, using the methods of this chapter and and considering 
time direction in which hops away from T. (and drift up) one can prove 
easily that contribution of Xreg.bound to the remainder is O^fj,"^^" h^^ at least 
as /i < h^-''. 

Consider zone Xdeg.inn- Then analysis of the previous chapter gives 
us remainder estimate 0(^fi~^^'^h~^ + {ix^^'^hY^^h'^y Recall that the last 
term is responsible for the correction which would come from near-periodic 
trajectories if 1/1/ = const. Now it is time to use the last non- degeneracy 
condition 



(5.17) 



\VtW\ 



where V* and W = ~V*\y: are introduced by formulae (14.2.28) and (14.2.29) 
[Ivr2]. Then using a canonical form of the previous chapter near T. (and 
thus ignoring presence of the boundary) and considering evolution of ^2 one 
can take T^ = h^~^i~^ with £ x |xi| (and this takes in account the presence 
oidX). 

As T* < e/i""^/*^ one can then ignore all periodicity, which proves that 
the contribution of zone X^^g H {xi > i} to the remainder is 0{ix~^^'^h~^); 
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meanwhile, the contribution of zone Xdegjnn H {xi < £} to the remainder 
is 0{i^-^/''h-^ + {fi^/^hy/^h-^i) with I = C maxifi^/^h^^^ fi-^/"). By our 
standard technique we can actually make 5 = in the estimate and then we 
arrive to the remainder estimate 0{^iJi~^/^h~^ + [ji^^'^hy^^h^^Y In particular, 
remainder estimate is 0{n^^^'^h^^) as /x < h~^^/^. 

5.4.3 Strong magnetic field case 

Here comes the difficult part: 

Problem 5.8. (i) In frames of condition (5.3) prove that if 

(5.18) 16 - r]lW^'^\ <p< t'^, with pl>h 

then one can take T^, x hl~^ and T* x 1 and then contribution of such 
element to the Tauberian remainder does not exceed 

(5.19) Cfi'^/'^pih-^ (1 + p^/" T,) T*'-^ X Cfi-^/''h-^pi{l + p^'^hr^) 

and therefore the total contribution to the Tauberian remainder of elements 
of this type with (i > max{h^/^, p'^/") is 0(/x-^/^/?~^).; 

(ii) Prove that if 

(5.20) \C2-vl^^^^\^P>^^^^. with pi >h 

then one can take T* x hp^^ and T* ^ p and then contribution of such 
element to the Tauberian remainder does not exceed 

(5.21) Cp-'/'pih-' (1 + /i^/" T,) T*-~' X Cp-'/'h-'piil + p'^''hp~^)p-' 

and therefore the total contribution to the Tauberian remainder of elements 
of this type with p > max{h^/^ , p-'^/^") is Oip-^^^h-^ + \ \ogp\); 

(iii) Provide the same estimate for contribution of the remaining small 
zone {p < max{h^/^, p~^/^'^), Cp~^/^ < ^ < h^^^,pi < h}; combining with 
estimate for Xcomer we arrive to the same estimate 0{p~^/^h~^ + /?~^); 

(iv) Improve this estimate to O^p'^'^h"^ + 1); 
Also looks challenging the next step: 
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Problem 5.9. Pass from asymptotics with the Tauberian principal part to a 
more exphcit ones. 

There should be usual Magnetic Weyl term, a term associated with the 
boundary and a term associated with degeneration albeit we do not expect 
any mixed boundary- degeneration term. 

5.4.4 Superstrong and hyperstrong magnetic field 
cases 

In the same way one needs to consider modified Schrodinger operator (with 
3/i/)F* subtracted where F"^ is "F but without absolute value") and taking 
into account remark 5.7 follow the same path: 

Problem 5.10. Derive spectral asymptotics for the modified Schrodinger 
operator: 

(i) Consider the case of superstrong magnetic field: eh~^ < /i < eh~^:, 

(ii) Consider the case of hyperstrong magnetic field: /i > eh~^ . 



6 Appendices: eigenvalues of L(r]) 
6. A Basic properties 

Recall that Xo^nij]) and A|\i,n(^) n = 0, 1, ... are eigenvalues of^"^ 

(6.1) L{r,) = -dl + {x + r^f 

on ]R~ with Dirichlet and Neumann boundary condition respectively at 0, 
Un{x) = Vn{x + Tj, Tj) are (real-valued) eigenf unctions such that \\vn\\ = 1. 

Proposition 6.1. AD,n(^) o-nd \N,n{v)> n = 0,1, 2, ... are real analytic func- 
tions with the following properties: 

(i) ^D.n{v) o,re monotone decreasing for rj G M; AD,fi(^) /^ +oo as rj ^ — oo; 
>^D.n{v) \ (2n + 1) asT]^ +cx); Ad,„(0) = {An + 3); 



20) 



Obviously, transformation x h^ — x reduces (1-26) to this one. 
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(ii) AN,n(^) are monotone decreasing for rj G M^; X^.n^n) Z' +00 asrj ^ —00; 
AN,.(r7) < (2n + 1) as r^ > (2n + 1)^; An,.(0) = (4n + 1); 

(m) XH.niff) < Xo.niv) < XN.{n+l){v)i ^D.niv) > (2"+!), ^N.niv) > (2"-!)+. 

Proof. Note that Xo.niv) ^^^ An „(?]) are obtained by variational procedure 
from quadratic forms 

(6.2) Q{u)=[ (\d,u\^ + {x + r]f\uf)clx, Qo{u)=[ \uf dx 

Jx<0 ^ ' Jx<0 

respectively with Dirichlet or no boundary conditions at x = 0. Alternatively 
we can consider 

(6.3) Q'{u)= f (\d^u\^ + x^\u\^)dx, Q^{u)= f \u\^ dx 

J x<r] J x<r} 

respectively with Dirichlet or no boundary conditions at x = rj and we will 
shift notations between these two. 

We know eigenvalues A^ = (2n + 1) and eigenfunctions of this operator on 
M (i.e. as ?7 = — oo). Further, we know that eigenvalues of /-d(0) and /-n(0) 
are (2m + 1) respectively with odd and even m we arrive to all assertions in 
(i), and to all assertions in (ii) except the third one. 

Consider Vk{x) which are eigenfunctions of harmonic oscillator on M, 
k = 0, ... , n. As u belongs to their span 

Q-{u) = Q{u) - Q+{u) < (2n + l)Qo{u) - Q+{u) = 

(2n + l)Qo{u) + (2n + l)Q^{u) - Q+{u) 

where Q, Qq denote the same forms on M and Q^, Qq on {x > rj}. 

Let rj > (2n + 1)2. Obviously (2n + 1)Qq{u) < Q^{u) and therefore 
Q~{u) < (2n + 1)Qq[u) for u belonging to some (n + l)-dimensional space; 
then Xu.niv) < (2n + l). 

In (iii) one needs only to prove that Xoniv) < AN,(n+i)(^)- However in 
(n + 2) dimensional span of V[^fi{x), ... , V[^(^n+i){x), functions vanishing at 
X = rj form (n + 1) dimensional subspace on which Q^{u) < Xf^ (^n+i)Qo {'^) 
(strict inequality is due to the fact that only on v = vt^_(^n+i){x) equality can 
be reached and this function does not belong to the subspace in question 
(as then v[ri) = dv[ri) = which is impossible). D 

So far we have not proved that AN,f,(^) — ?■ (2n + 1) as 77 — )■ +00; it will 
be done later. 
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6.B More properties 

We follow M. Dauge-B. Helffer [DH] here. 
Consider equation on x < 

(6.4) {L{r^)-\,{r^))u, = Q. 
Differentiating by t] and multiplying by u^ we arrive to 

(6.5) dr^\n{'n) = {{d^L{r]))Un, Un) = 2((x + r])Ur„ Un) = 

([5x, (x + r]f]un. u,) = 2((x + r]fun, d,U„) + V^\Un{0)f = 

{r]^-\n{v))\Un{0)f-\dMO)f- 

In particular, 

(6.6)d drjXo.niv) = -\dxVD.n{v)f 

and 

(6-6)n d^XN.niv) = [if - ^N,f,(^))|^N,n(^)|^- 

We refer to them as Dauge-Helffer formulae [DH]. 
We conclude immediately that 

Proposition 6.2. (i) driXo.niv) < 0/ 
(a) drfXti^niv) ^ ^/ (^''T'd only if 

(6.7) X^,niv) I V^. 
Differentiating (6.6),^ we arrive to 

(6.8) d^^XN,n{v) = 

Therefore as 9^AN,n(^) = in virtue of proposition 6.1 

(6.9) d^^XN,n{v) = HvN.n{v)\^ 

and therefore it is non-degenerate minimum. So, 
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Proposition 6.3. (i) AN,n(^) has a single stationary point^^^ r]n, it is non- 
degenerate minimum, and in this point (6.7), (6.9) hold. 

(a) In particular, {Xu,n{v) ~ 1^) has the same sign as [rjn — rj). 
Therefore X^.n+iiVn) - vl > X^.tiiVn) - vl = ^ implies tliat 

(6.10) r]n+l > Tjn. 

6.C Estimates of X*,n{v) as 77 — > +00 

So we must investigate properties of Xniv) as r^ — )■ +00. We know tliat tlien 
XD,n{v) > (2n + 1) > Xu^niv) and Xo.niv) tends to (2n + 1). 
Consider Vj{x). Note that 

r+ca 

(6.11) / vj{x)vkix)dx=0{rj'+^-^e-'^') 

(non-uniformly with respect to 7, k) and a,s J = k this is an exact magnitude. 
Therefore 

(6.12) {vj,v,) = 8j, + 0{r/+'-'e-'>") 

as our space is ^^((—00,77)). 
Consider bihnear form 

{u'v' + x^uv) dx 

■00 

(where ' means dx)- One can rewrite it as 

/v 
[-u" + x'^u)vdx+ u'{r])v{r)). 
-00 

Then 

Q~{vj, Vk) = (2j + l){vj, Vk) + Vj{r])vk{r]). 



'^^' And it must have one due to proposition 6.1. 
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Thus if we consider u = J2k<n '^k^k with «/< G M such that J2k<n <^| ^ 1 
then 

Q-{u) - (2n + l)Q^{u) = J2^{k- n)al + 0{e-^\"^'^') 

k<n 

and this expression is negative provided a^ > e~^ ^^n+k+i f^^. gQj^g k < n; 
otherwise we can assume that «„ = 1 and 

Q-{u) - (2n + l)Q^{u) ^ 5]((?~(t)/c) - (2n + l)Q,{vk))al < 

k<n 

{Q-{vn)-{2n + l)Q^{v,))=v'Myn{v) mod 0(e-|^') 

and we know that v'^{v)^ n{v) is negative and of magnitude 7^2n+ig-r? g^g 
T) > Cn and therefore due to variational principle 

XnAv) < (2n + 1) - eT^^^+^e-"' as r^ > c„. 

Meanwhile let us replace Vk{x) by Vk{x) = Vn{x) + (f)k, <Pk = ~Vk{rj)e'^^~'^ 
satisfying Dirichlet conditions a,s x = r]. All the above conclusions except 
calculations of Q{vn) remain true for these functions as well and for them 

Q~{v,)-{2n + l)Q,{v,) = 

v'n{v)vn{v) + 2vM(l),{r]) + (?-(0„) - (2n + l)(?o"(0n) = 

- vMvniv) + vliv) {V + 0(r^)) X r^2"+ie-^= 
and therefore 

^dAv) < (2n + 1) + CT^^^+^e-"' as r/ > c„. 

On the other hand, consider ^ 3 u supported in (— oo, rj] and orthogonal 
to -uo, ■■■ ,Vn-i. Then (modulo constant factor) u = Vn{x) + v where v is 
orthogonal to Vn in M as well. 

Note that 

Q-{u) = Q{u) = Q{vn) + Q{v) = (2n + 1) + Q{v), 

and Qo{u) = 1 + ||v||^, and 

Q-{u) - (2n + l)Qo{u) = Q{v) - (2n + l)Qo{v) > 2Qo{v) > 2Q+{v) = 

2Q+{v„)>er]^"+'e-^\ 
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So, any (n + l)-dimensional subspace of ^^((— oo, t])) contains an element 
with (?~(u) > ((2n + 1) + erf"^^e^'^ )Qq{u) and we arrive to the left 
inequality of 

(6.13) enr]^"+^e-^' < XoM " (2n + 1) < Q^^'+'e-"' as r] > c, 

as the right one has been already proven. 

Finally, estimate An „ from below. Consider Vn{x) and extend it as 
Vn{Ti)e2^''' ~^ ) to [i], oo). Let u be a resulting function. Note 

Q+{u) - (2n + l)Q+{u) = \v„{v)\'{v - 2n<P{r])). 
with ^(77) = /°^ e^ "^ c/x and then as n = 

An,„(77) - (2n + 1) = (?-(i;,) - (2n + l)(?o"K) = 

(?(ty) - (2n + l)(?o(tv) - (?+(u) + (2n + l)(?o+(t^) > 

- K{v)\'{v - 2ncD(r7)) = -^(^ - 2n<i>{r])){r]' - XM.nr'd,\M.n{v) 
due to ( 6. 6 )|vg. Therefore 

(2n + 1 - AN,.(r/))"'9,AN,n(r/) > 2(77^ - AN,.(r/)) (r/ - cl>(r/))"\ 

Note that 0(77) = ^r^^^ + 0(77^^) and An < 2n + 1. Then the right-hand 
side of the last inequality is greater than 2{ri — [n + l)ri'^) — 0{ri~^) and 
integrating we conclude that 

log(2n + 1 - XnAv))'^ >rf- 2(n + 1) log 77 - log C 

and 

2n + 1 - An,.(77) < a„(r7) := C,e-'^" r,^"+^ 

which is just one factor t] greater than we got as an lower estimate for 
(1-A„). 

Finally, as AN,n+i(^) > AN,n+i(^n+i) > ^^oAnn+i) > (2n + 1) we conclude 
that \N,n+i{v) ^^^ ^N,n{v) ^'^'^ disjoint as ?] > Q and therefore as Q~{vn) < 

Q~{Vn) + CnCTn WC COUcludc that Q^{Vn — Vn) < CnCrnlr]): 

(6.14) \\d>,{Vn - Vn)\\ + \\x{Vn - Vn)\\ < Cn(Jn['nY 
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where norms are calculated in =Sf^((— oo.r])) and then 

(6.15) \{vn-Vn)\<C,an{r]Y^- 

In particular jt^nl^)! < Qo"n(^)2- 

So far it was proven only for n = 0. However we apply induction: as it 
was proven ioi k < n with fast decaying (Tk[f])- Then {vn,Vk) = Oiak^r])^^ 
ioi k < n and extending Vn as before we arrive to (tv, Vk) = 0[an-i{r])2^ for 
all k < n. 

Then Q{u) — (2n + 1)Qq{u) > —Can-i and we arrive to 

An,. - (2n + 1) > -^(77 - 2n<i'{'n)) {rf - \^,n{v)y'dr,\uAv) " Cr]^"e-^" 

which implies the same conclusions. 

We also conclude that (6.14) and (6.15) hold for Vn replaced by vo,n and 

that \d:^VD,n{r])\ < CnTjaniv)^- 

6.D Asymptotics of X*,n{v) as 77 — > +oc 

Here will reproduce results of section 5 of C. Bolley-B. Helffer [BH] and 
generalize them to arbitrary j and to Dirichlet problem (definitely their 
methods could do the same). 

In contrast to this paper we will not build quasi-modes but use Dauge- 
Helffer formulae. Let 

(6.16) e,,n{v) = HKn{v)-^n-l) 
with * = D, N respectively. Then from 

L{v)v*.n{v) - (2n + 1 ± e*,„(r7))t;*,„ = 0, L{r])v„ - (2n + l)t}„ = 

we conclude that 

Te*,n{v){'^*.n.Vn) = v'„{r])v ,{r]) - v,{r])v'^{r]) 
and as modulo 0(e^2^' r]'^^ {v*,n:Vn) = 1 then 
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modulo 0(e^2'' r]^^. Then modulo sign 

1 

1 

and modulo 0(^6 2^ r]'^) 

and 

(6.17)n e^'niv) = fiv^ - 2n - l)(i):(r7))-' d^ 

and finally we arrive to 
Theorem 6.4. As rj ^ +oo 

(6.18)d AD,n(r7) = 2n + 1 + [j\v,{r^)y^ dr^y', 

(6.18)^ An,.(77) = 2n + 1 - (y"'(772 _ 2„ _ l){vMy' dri)'' 

modulo 0[e~2^ r^'^) with K = K{n). 

Remark 6.5. One can see easily from (6.18)p, (6.18)|vj that 

k>0 

with the leading coefficient c'^ no = |/t^^ where Kn = (n!)^27r^42"/^ coefficient 
at x" at VnS'^ /^ and thus 

(6.19) e~^'eUv) =~ Yl ^*.nW"^'-^' 

k>0 

where 

on+l 

2 ^ 



(6.20) Q,„,o = 2< 



nivvr 
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Now we can estimate recurrently higher-order derivatives of e*,n(^), 
v^,n{x, 1]) as X = 1] arriving to 

Corollary 6.6. 

(6.21) \&X,,n{v)\^V^"^^^'e~''' 

as 1] > CnJ. 

6.E Remarks 

Remark 6. 7. While we do not need it one can prove easily that 

(6.22) A,,„(r^)~r^2^5^Q,„,,7^i(^-2.) 

as r] — 7- — oo with c*,n^o eigenvalues of operator Mirf) := (D^ + 2x|r7|) on IR+ 
with the corresponding boundary condition (D or N) as x = and thus are 
defined through zeros of Airy function or its derivative. 

Remark 6.8. (i) For operator D^ +x^ on {x < 7]} albeit with the boundary 
condition [u' + au){vi) = with < a < oo one can see easily that (6.6)^, 
(6.17)pj are replaced by 

(6.6)^ d^\n{r], a) = {rj^ - c? - \n{r], a)) \vn{r], a)|^ 

(6.17)^ e-\r], a) ^ j\rf - 2n - 1 - a'){v'^[r^, a))-^ dr^ 

(6.23) 

respectively. 

(ii) One can prove easily that \n{i]', ci) monotonically depends on a and 

Xn{rj; 0) = Xu.niv), ^n{V: ") /" ^D.niv) aS tt -> +00. 

However, in virtue of (i) for all a > there is a single extremum which 
is a nondegenerate minimum at rjn = rin{a) > such that rj^ = \n{f]n) + o;^ 
and as q; — )■ +oo therefore rin{a) ~ a^ + (2n + 1) with a small negative error 
(see remark 6.8). 

Therefore as a — > +oo these minima (see blue line on figure 12 for a = 0) 
become more shallow and move to the right and blue-like lines here can 
approximate red-like lines exactly as on the graphs of e"'' — [rj^ + a^)~-^.). 
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Figure 11: Schematics of ({rf) = Xnii]; a) — (2n + 1). 



Conjecture 6.9. For each n function Xo.niv) ^-^ convex and d^Xo^niv) > 
while \N,n{v) has a single non- degenerate inflection point tj'^ and d^Xoniv) 



^ 



as rj ^ rj'^. 
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Dirichlet versus Neumann Boundary Conditions 
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Figure 12: Plots of AN,f,(^) and Xo.niv)- Calculated by Matlab: Courtesy 
of Dr. Marina Chugunova. Here rj is replaced by rj 
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